D  D  C 

-  '"-mnnE 

APR  1 2  1967 


iU 

C 


PRINCETON  UNIVERSITY 

DEPARTMENT  OF 

AEROSPACE  AND  MECHANICAL  SCIENCES 


mm  ©« 


itALimo  KfcOUfciTCRS  MAT  0UTA1N  AUOITIONAL  COPlfcS  MOM 

rm  r.Tft.?r  rocuww'w  center,  all  others  should 
'  \f  -iHOUSE  FCR  FEDERAL  SCIENTIFIC  AND 


Best 

Available 

Copy 


tap 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
AFOSR  REPORT  NO.  66-2725 

i 


EXPERIMENTAL  INVESTIGATION  OF  HIGH-FREQUENCY 
LONGITUDINAL  COMBUSTION  INSTABILITY  IN 
GASEOUS  PROPELLANT  ROCKET  MOTORS 

BY 

CRAIG  THOMAS  BOWMAN 

AFOSR  Contract  AF49  (638)1268 
Project-Task  9711-01 

January  1967 

K 


Approved  by: 

Irvin  Glassman 
Professor  of  Aeronautical 
Engineering 

Principal  Investigator 


Qualified  requestors  may  obtain  additional  copies 
from  the  Defense  Documentation  Center.  All  others 
should  apply  to  the  Clearinghouse  for  Federal 
Scientific  and  Technical  Information. 

Reproduction,  translation,  publication,  use  and 
disposal  in  whole  or  in  part  by  or  for  the  United 
States  Government  is  permitted. 


Guggenheim  Laboratories  for  the  Aerospace  Propulsion  Sciences 
Department  of  Aeroepace  and  Mechanical  Sciences 

Technical  Report  784 

Princeton  University 
Princeton,  New  Jersey 


A 


A  dissertation  submitted  to  the  Department  of 
Aerospace  and  Mechanical  Sciences  of  Princeton 
University  in  partial  fulfillment  of  require¬ 
ments  for  the  degree  of  Doctor  of  Philosophy. 


Guggenheim  Laboratories  for  the  Aerospace  Propulsion  Sciences 
Department  of  Aerospace  and  Mechanical  Sciences 
Princeton  University 
Princeton,  New  Jersey 
1966 


V 


ii 

ABSTRACT 

The  properties  of  the  longitudinal  mode  of  high-frequency  com¬ 
bustion  instability  in  a  rocket  motor  burning  premixed  gaseous  propellants 
have  been  determined  experimentally.  The  experimental  observations  have 
been  compared  with  the  results  of  a  non-linear  instability  theory  (23)  " 
based  on  a  chemical  kinetic  driving  mechanism.  It  has  been  shown  that  the 
theory  gives  the  correct  waveform  for  longitudinal  instability.  In  addi¬ 
tion  the  theory  gives  the  correct  qualitative  dependence  of  the  stability 
limits  on  the  mean  combustion  temperature  and  activation  energy.  The 
theory  also  gives  the  correct  qualitative  dependence  of  the  instability 
strength  on  the  mean  combustion  temperature  and  mean  Mach  number  of  the 
combustion  gases.  Harmonic  mode  longitudinal  instabilities  have  been 
observed,  and  a  criterion  for  their  appearance,  consistent  with  a  chemical 
kinetic  driving  mechanism,  has  been  suggested.  The  dependence  of  the 
stability  limits  and  instability  strength  on  the  mean  combustion  pressure, 
combustion  chamber  lei  gth,  injector  and  exhaust  nozzle  have  been  deter¬ 
mined  experimentally.  These  effects  are  not  included  in  the  instability 
theory.  Plausible  explanations  for  the  above  effects,  consistent  with  the 
theoretical  model,  have  been  suggested.  Based  on  the  general  agreement 
between  the  experimental  observations  and  theoretical  results  it  is  con¬ 
cluded  that  the  appropriate  driving  mechanism  for  high-frequency  longi¬ 
tudinal  gas-phase  combustion  instability  is  chemical  kinetics. 
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CHAPTER  1.  INTRODUCTION 

The  phenomenon  of  combustion  instability  in  rocket  motors  is 
characterized  by  periodic  oscillations  in  the  combustion  pressure. 

These  pressure  oscillations  generally  result  in  erratic  operation,  and 
frequently  cause  a  catastrophic  failure  of  the  rocket  motor.  It  is 
desirable,  therefore,  to  understand  the  nature  of  the  instability 
phenomenon  and  to  derive  design  criteria  for  preventing  unstable 
combustion. 

In  general  three  basically  different  types  of  combustion 
instability  are  recognized.  The  first  of  these  is  termed  "low- frequency 
instability"  or  "chugging",  and  i6  characterized  by  oscillation  fre¬ 
quencies  in  the  range  of  20  to  200  cps.  Low- frequency  instability 
results  from  an  interaction  between  the  combustion  process  and  the  pro¬ 
pellant  feed  system.  Studies  by  various  investigators  (l-12)  have 
resulted  in  a  satisfactory  understanding  of  low-frequency  instability, 
and  have  provided  reasonably  accurate  guidelines  for  preventing  its 
occurence. 


Ibe  second  type  of  instability  is  termed  "  interned  late  - 


♦Numbers  in  parentheses  indicate  references  listed  in  the 


references  section. 


frequency  instability",  and  is  characterized  by  oscillation  frequencies 
in  the  range  of  200  to  600  cps.  Intermediate-frequency  instability  results 
from  the  production  of  entropy  vaves  in  the  combustion  chamber  and  the 
interaction  of  these  waves  with  the  exhaust  nozzle.  Several  studies  of 
the  intermediate-frequency  instability  phenomenon  have  been  made  (10,13). 
Intermediate-frequency  instability  is  seldom  obser  ed  in  practice. 

The  third  type  of  instability  is  termed  "high-frequency 
instability"  or  "screaming",  and  is  characterized  by  oscillation 
frequencies  generally  in  excess  of  600  cps.  High-frequency  instability 
results  from  an  interaction  between  the  combustion  process  and  the 
acoustic  properties  of  the  combustion  chamber,  and  exists  in  several 
forms  which  depend  on  the  acoustic  mode  of  the  chamber  excited  by  the 
combustion  process.  Thus  a  distinction  is  made  between  longitudinal  and 
transverse  mode  instabilities.  High-frequency  combustion  instabiltiy  has 
received  a  great  deal  of  attention  in  recent  years.  Experimental  studies 
by  various  investigators  (9-12»  14-17)  have  helped  to  define  the  nature 
of  the  phenomenon.  Theoretical  studies  based  on  various  simplified 
models  have  been  made  by  Crocco  (5,  8,  17-19) >  Priem  and  Guentert  (20), 
Culick  (21),  Hart  and  McClure  (22)  and  others.  These  studies  have 
provided  a  partial  understanding  of  the  phenomenon  together  with 
suggestions  on  how  to  improve  the  stability  characteristics  of  a  rocket 
motor.  The  fact  that  a  complete  theory  of  high-frequency  combustion 
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lnstabillty  in  solid  and  liquid  propellant  motors  is  not  yet  available 
is  largely  due  to  uncertainties  in  the  details  of  the  unsteady  multi¬ 
phase  burning  process. 

From  a  theoretical  standpoint  it  is  desirable  to  separate  the 
complex  multi-phase  burning  process  into  a  series  of  simpler  processes 
(eg.  atomization,  vaporization,  mixing,  reaction).  The  contribution 
that  each  individual  process  makes  to  the  overall  instability  phenomenon 
then  can  letermined.  If  it  were  possible  to  couple  the  contributions 
of  the  individual  processes  an  overall  instability  theory  could  be  ob¬ 
tained.  Such  a  coupling  may  not  be  possible  but,  the  investigation  of 
the  individual  processes  may  provide  certain  stability  parameters  which 
would  not  be  available  from  a  simplified  combustion  model. 

The  present  work  is  directed  toward  understanding  the  contri¬ 
bution  of  the  gas  phase  to  high-frequency  combustion  instability.  The 
investigation  employs  a  rocket  using  premixed  gaseous  propellants.  This 
"gas  rocket"  eliminates  those  phenomena  peculiar  to  solid  and  liquid 
prcoellant  motors,  and  any  instability  encountered  can  be  attributed  to 
phenomena  in  the  gas  phase.  By  its  very  simplicity  the  gas  rocket 
readily  lends  itself  as  a  model  for  the  theoretical'  treatment  of  the 
problem. 


Simplicity  is  only  one  reason  for  investigating  gas-phase 
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combustion  Instability.  Experimental  studies  (l4,  17,  24)  have  shown 
that  there  are  certain  similarities  between  high-frequency  instability 
in  gas  rockets  and  high-frequency  instability  in  solid  and  liquid 
propellant  motors  (eg.  pressure  waveform).  In  addition  it  is  possible 
that  the  gas  phase  plays  an  important  role  in  the  Initiation  of  high- 
frequency  oscillations  in  solid  and  liquid  propellant  motors.  Hence 
an  understanding  of  gas-phase  instability  may  provide  some  insight  into 
the  more  complicated  instability  phenomena  observed  in  solid  and  liquid 
propellant  motors. 

In  this  thesis  the  experimentally  observed  properties  of  the 
longitudinal  mode  of  high-frequency  combustion  instability  in  a  gas 
rocket  are  presented.  A  comparison  of  the  experimental  observations 
with  the  results  of  an  instability  theory  (23)  based  on  a  model  closely 
approximating  conditions  in  a  gas  rocket  is  made.  From  this  comparison 
some  conclusions  concerning  the  driving  mechanism  for  high-frequency 
Instability  in  gas  rockets  are  drawn,  and  some  comments  on  the  role  of 
the  gas  phase  in  high-frequency  instability  in  solid  and  liquid  pro¬ 
pellant  motors  are  made. 
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CHAPTER  2.  EXPERIMENTAL  APPARATUS 

A.  Gas  Rocket  System 

From  experimental  studies  of  high-frequency  combustion 
instability  in  liquid  propellant  rocket  motors  (10-12,  14-17)  it  is 
known  that  the  longitudinal  mode  is  strongly  dependent  on  the  length 
of  the  combustion  chamber,  the  mixture  ratio,  the  propellant  combination 
and  also  to  a  lesser  extent  on  the  injector  and  the  nozzle.  The  gas 
rocket  system  was  designed  so  that  these  five  parameters  could  be 
changed  easily.  Two  basic  sub-systems  make  up  the  gas  rocket  system- 
the  propellant  feed  system  and  the  rocket  motor  (consisting  of  a  mixing 
chamber,  injector,  combustion  chamber  and  exhaust  nozzle). 

1.  Propellant  Feed  System 

The  purpose  of  the  propellant  feed  system  (Figure  1)  is  to 
supply  the  desired  propellant  flow  rates  to  the  rocket  motor.  The  feed 
system  is  capable  of  simultaneously  supplying  as  many  as  three  different 
gases  to  the  motor  (The  need  for  this  capability  will  be  discussed  in 
Chapter  4).  Flow  metering  is  accomplished  using  critical  flow  orifices 
of  the  type  shown  in  Figure  2.  A  total  of  ten  orifices  with  different 
throat  diameters  are  available  so  that  a  wide  range  of  flow  rates  can 
be  obtained.  The  critical  flow  orifices  were  calibrated  using  a  water 
displacement  technique  (Appendix  A). 
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2 .  Rocket  Motor 

The  rocket  motor  sub- system  consists  of  a  mixing  chamber,  an 
injector,  a  combustion  chamber  and  an  exhaust  nozzle.  A  schematic  of 
the  rocket  motor  is  shown  in  Figure  3*  and  a  photograph  of  the  motor 
mounted  on  the  test  stand  16  shown  in  Figure  U.  The  individual 
components  of  the  rocket  motor  are  described  below. 

Mixing  Chamber 

The  mixing  chamber  (Figure  5)  provides  for  the  mixing  of  the 
gaseous  propellants  prior  to  their  injection  into  the  combustion 
chamber.  Primary  mixing  of  the  fuel  and  oxidizer  is  provided  by  the 
impingement  injector.  Provision  is  made  for  introducing  a  third  gas 
Just  downstream  from  the  impingement  injector  through  a  radial- type 
showerhead  injector.  Final  mixing  is  obtained  by  passing  the  crudely 
mixed  gases  through  a  tightly  packed  bed  of  stainless  steel  balls. 

Experiments  performed  by  Zucrow  (2k)  using  premixed  and  un¬ 
mixed  gaseous  propellants  show  that  the  degree  of  mixing  can  influence 
the  instability  phenomenon.  It  is  of  interest,  therefore,  to  consider 
the  extent  to  which  the  propellant  gases  are  mixed  by  the  mixing  chamber. 
To  thi'.  end  an  experiment  was  carried  out  in  which  given  flow  rates  of 
an  oxidizer  (0^)  and  a  fuel  (H^)  were  supplied  to  the  mixing  chamber. 


Samples  of  the  gas  mixture  at  a  station  Just  downstream  from  the  mixing 
chamber  were  obtained.  These  samples  were  analyzed  using  a  gas 
chromatograph  to  determine  the  relative  amounts  of  Og  and  Hg  (le. 
mixture  ratio).  This  mixture  ratio  was  then  compared  to  the  input 
mixture  ratio  (determined  from  the  known  flow  rates  of  Og  and  Hg)  to 
obtain  an  Indication  of  the  degree  of  mixing.  In  the  experiment  three 
different  mixture  ratios  were  used;  five  samples  were  taken  for  each 
mixture  ratio.  The  results  of  this  experiment  are  shown  in  Figure  6. 
The  scatter  in  the  data,  as  Indicated  by  the  vertical  bars,  is 
approximately  *  This  Is  of  the  same  order  as  the  uncertainty  In 
the  known  flow  rates  (Appendix  B).  It  Is  concluded,  therefore,  that 
the  mixing  chamber  provides  for  complete  mixing  of  the  gaseous 
propellants. 

Since  the  propellants  are  premixed  before  Injection  Into  the 
combustion  chamber  there  exists  the  possibility  that  combustion  (which 
normally  occurs  In  the  combustion  chamber)  can  propagate  upstream 
through  the  injector  and  occur  In  the  mixing  chamber  itself.  The 
tendency  for  such  a  "flash-back"  of  combustion  cam  oe  considerably 
reduced  by  a  correct  choice  of  injector  and  starting  procedure. 
Unfortunately  even  with  these  precautions  flash- back  does  occur  on 
occasion  (partic»Uarly  when  operating  with  a  propellant  combination 
for  the  first  time).  To  prevent  damage  to  the  mixing  chamber  by 


flash-back  two  safety  devices  were  Incorporated  in  the  chamber  design. 
Hie  first  device  is  a  burst  diaphragm,  designed  to  relieve  the  over¬ 
pressure  associated  with  combustion  in  the  mixing  chamber.  Hie  second 
device  is  an  electro- me  chani  xl  system  (Figure  7)  which  automatically 
shuts  off  the  oxidizer  flow  when  flash-back  occurs  (25) .  Hie  system 
consists  of  a  sensing  circuit  which  actuates  a  solenoid  valve  which  in 
turn  controls  a  pressure-actuated  valve  in  the  oxidizer  line  upstream 
from  the  mixing  chamber.  Hie  sensing  elements  are  photoelectric  cells 
mounted  in  the  wall  of  the  mixing  chamber,  and  sense  the  combustion 
generated  light  through  plexiglas  windows.  After  a  shutdown  due  to 
flash-back  fuel  continues  to  flow  to  prevent  any  remaining  oxidizer 
from  reacting  with  the  metal  components  of  the  mixing  chamber. 

Hie  mixing  chamber  is  also  fitted  with  a  pressure  tap  so 
that  the  static  pressure  of  the  gases  in  the  chamber  can  be  monitored 
continuously. 

Injectors 

After  being  thoroughly  mixed  in  the  mixing  chamber,  the 
gaseous  propellants  flow  through  an  injector  into  the  combustion 
chamber.  Hie  injector  is  essentially  a  circular  metal  plate,  with 
some  arrangement  of  injection  orifices,  fitted  in  the  downstream  end 
of  the  mixing  chamber  (Figures  3  and  5)»  To  determine  the  effect  of 
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the  injection  pattern  on  the  instability  phenomenon  experiments  were 
carried  out  using  three  different  types  of  injectors.  The  first  type 
of  injector  is  a  "porous"  plug  composed  of  sintered  stainless  steel 
powder  with  a  mean  pore  size  of  five  microns  (Figure  8) .  The  porous  plug 
injector  provides  nearly  one  -  dime  nsional  injection  of  the  propellants 
into  the  combustion  chamber  and  hence  should  provide  experimental  data 
which  can  be  correlated  with  one-dimensional  instability  theories. 

To  determine  the  effect  of  deviations  from  one -dimensional 
injection  on  the  instability  phenomenon  two  types  of  injectors  which 
produce  combustion  z-»ne  recirculation  were  used.  The  first  type,  termed 
a  'fehowerhead"  injector,  is  a  copper  plug  with  injection  orifices  parallel 
to  the  longitudinal  axis  of  the  combustion  chamber.  Showerhead  injectors 
characteristically  produce  weak  recirculation  patterns;  hence,  slight 
deviations  from  one-dimensional  injection  are  to  be  expected.  Three 
different  showerhead  injectors  were  used  in  the  instability  studies-  10 
holes,  31  holes  and  49  holes  (Figures  8  and  9)»  The  second  type  of 
"recirculating"  injector  is  a  copper  plug  with  two  sets  of  impinging 
triplets  (Figure  10) .  Impinging  injectors  produce  strong  recirculation 
patterns;  hence,  fairly  significant  deviations  from  one-dimensional 
injection  are  to  be  expected.  The  total  flow  areas  for  the  three 
showerhead  injectors  and  the  impinging  injector  were  approximately  equa.., 
and  were  chosen  so  that  the  injection  flow  was  choked  under  all  conditions 
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of  operation.  Therefore  any  instability  encountered  using  these  injectors 
cannot  be  attributed  to  variations  in  the  propellant  injection  rate.  It 
vas  not  possible  to  obtain  choked  injection  flow  using  the  porous  plug 
injector. 


Combustion  Chamber 

The  cylindrical  combustion  chamber  has  an  internal  diameter  of 
3.8  cm.  This  diameter  was  chosen  so  as  to  discourage  the  occurence  of 
the  transverse  modes  of  combustion  instability  (As  will  be  seen  in 
Chapter  ^  transverse  modes  do  occur  on  occasion  near  stability  limits). 

The  length  of  the  combustion  chamber  can  be  varied  from  10  cm  to  200  cm 
by  using  various  combinations  of  nine  12.5  cm  sections  and  two  50  cm 
sections  together  with  one  of  three  n  zzle  assemblies  (described  in  the 
following  section).  A  typical  combustion  chamber  section  is  3hown  in 
Figure  11,  The  design  of  the  chamber  sections  took  into  consideration 
the  large  heat  transfer  rates  associated  with  unstable  combustion  and  the 
requirement  of  extended  run  times  (30  seconds  or  more).  The  sections  have 
thick  copper  walls  which  serve  as  a  heat  sink  and  stainless  steel  liners 
to  reduce  heat  transfer  to  the  walls.  In  addition  the  exterior  wall  of 
the  combustion  chamber  was  spray-cooled  using  water. 

The  first  chamber  section  (ie.  the  one  Just  downstream  from  the 
mixing  chamber)  was  modified  to  accept  a  high-frequency  pressure  trans- 
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ducer,  a  spark  plug  (used  for  igniting  the  propellants),  a  steady-state 
pressure  transducer  and  five  thermocouples  (used  tc  obtain  temperature 
distributions  in  the  combustion  zone).  Two  other  chamber  sections  were 
each  modified  to  accept  a  single  high-frequency  pressure  transducer; 
hence,  data  from  as  many  as  three  high-frequency  pressure  transducers 
can  be  obtained  during  a  given  run. 

Exhaust  Nozzles 

Theoretical  studies  of  the  longitudinal  mode  of  high-frequency 
combustion  instability  (8,  21,  23)  indicate  that  the  instability  phe¬ 
nomenon  depends  on  the  mean  flow  velocity  of  the  combustion  gases  and 
on  the  geometry  of  the  exhaust  nozzle  (eg.  the  length  of  the  convergent 
section).  To  investigate  these  effects  experimentally  two  different  sets 
of  exhaust  nozzles  were  designed.  In  experiments  to  determine  the  effect  of 
the  mean  flow  velocity  on  the  instability  phenomenon  two  water-cooled  con¬ 
verging-diverging  nozzles  were  used  (Figure  12).  These  nozzles  have  iden¬ 
tical  convergence  angles  (60°)  and  different  contraction  ratio.}-  the  con¬ 
traction  ratio  (Cr  *  Achamber/Athroat)  for  one  nozzle  is  approximately  23 
and  for  the  other  nozzle  approximately  46.  In  the  experiments  the  exhaust 
flow  from  the  nozzles  was  choked  under  all  conditions  of  operation;  hence, 
the  mean  flow  velocities  for  the  two  nozzles  varied  inversely  with  the  con¬ 


traction  ratios.  To  study  the  effect  of  the  length  of  the  convergent  section  of 
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the  exhaust  nozzle  on  the  instability  phenomenon  experiments  were 
carried  out  using  the  smaller  converging-diverging  nozzle  and  an  un¬ 
cooled  ’ tungsten  "plug”  nozzle  (Figure  13).  These  tvo  nozzles  have 
approximately  the  same  flow  area  ,  but  have  converging  sections  of 
different  length  (The  plug  nozzle  approximates  the  "zero-length" 
nozzles  employed  in  some  theoretical  instability  models-  see,  for 
example,  Reference  23). 

B.  Instrumentation 

The  instrumentation  used  with  the  gas  rocket  Bystem  can  be 
divided  roughly  into  tvo  groups-  feed  system  instrumentation  and 
combustion  chamber  instrumentation. 

1.  Feed  System  Instrumentation  (Figure  l) 

The  desired  propellant  flow  rates  were  obtained  using  calibrated 
critical  flow  orifices.  For  sucn  orifices  the  flow  rate  is  determined  by 
the  stagnation  pressure  and  temperature  upstream  from  the  orifice.  Stag¬ 
nation  pressures  were  measured  using  steady-state  pressure  transducers 
(Statham  PG  731);  stagnation  temperatures  were  measured  using  copper- 
copper*  constantan  thermocouples.  The  output  signals  from  the  pressure 
transducers  and  thermocouples  were  recorded  on  Leeds  and  Northrup  record¬ 
ing  potentiometers.  To  assure  that  the  orifice  flov  was  choked  during  a 


given  run  the  pressure  downstream  from  the  orifice  was  monitored.  The 


pressure  in  the  mixing  chamber  was  also  monitored  to  determine  if  the 
injector  flow  was  choked. 

2.  Combustion  Chamber  Instrumentation 

Steady-State  Instrumentation 

The  steady- 6tai.e  or  mean  combustion  pressure  was  measured  by 
a  Flader  PS HD  pressure  transducer.  Mean  gas  temperatures  in  the  combus¬ 
tion  zone  were  measured  using  five  uncoated  platinum- plat inum+ 10$  rhodium 
thermocouples  with  platinum-6^  rhodium  shields  installed  in  the  first 
combustion  chamber  section.  Thermocouples  of  this  type  were  required 
because  of  the  high- temperature  oxidizing  environment  which  existed 
under  certain  conditions  of  operation  of  the  gas  rocket.  Because  of  this 
environment  the  thermocouples  were  not  equipped  with  radiation  shields. 

The  output  signals  from  the  five  thermocouples  and  the  Flader  transducer 
were  recorded  on  Leeds  and  Northrup  recording  potentiometers. 

Transient  Instrumentation 

The  amplitude  and  frequency  of  the  pressure  oscillations 
associated  with  unstable  combustion  were  measured  by  one  or  more  water- 
cooled,  strain  gage  type,  high-frequency  pressure  transducers  (Dynisco 
PT  49)  mounted  in  the  wall  of  the  combustion  chamber  so  that  the  pressure- 
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sensitive  diaphragms  were  flush  with  the  interior  wall.  The  output  sig¬ 
nals  from  the  Dynisco  transducers  were  recorded  on  a  seven-track  magnetic 
tape  recorder  (Ampex  S-356l).  For  the  purposes  of  data  reduction  the 
recorded  signals  were  displayed  on  a  visicorder  (Minneapolis-Honeywell 
1106).  The  only  filtering  of  the  transducer  output  is  that  inherent  in 
the  transient  pressure  recording  system.  The  effect  of  this  inherent 
filtering  on  the  transient  pressure  data  is  discussed  in  Appendix  B. 


CHAPTER  3.  THEORETICAL  CONSIDERATIONS 


A.  Historical  Background 

One  of  the  initial  works  to  appear  on  high-frequency  combustion 
instability  in  gas  rockets  was  due  to  Zucrow  and  Osborn  ( 31)  -  They  re¬ 
ported  the  results  of  an  experimental  study  of  the  longitudinal  mode  of 
combustion  instability  and  suggested  a  driving  mechanism  based  on  a 
coupling  of  the  pressure  oscillations  with  the  chemical  reaction  rate. 

In  this  driving  mechanism  a  small  localized  pressure  disturbance  (perhaps 
due  to  turbulence)  propagates  longitudinally  in  the  combustion  chamber. 

When  tills  pressure  disturbance  traverses  the  reaction  zone  it  increases  the 
burning  rate  (due  to  the  increased  pressure  and  temperature  associated  with 
the  disturbance).  The  increased  burning  rate,  in  turn,  reinforces  the 
disturbance.  If  the  energy  supplied  to  the  disturbance  (by  the  increased 
burning  rate)  is  greater  than  the  energy  lost  due  to  dissipation  and 
reflection  from  the  nozzle,  the  disturbance  will  be  amplified  during  each 
succeeding  oscillatior  until  dissipative  phenomena  within  the  wave  limit 
further  amplification.  If  the  energy  supplied  by  combustion  is  less  than 
the  energy  lost  the  disturbance  will  decay  in  time.  To  account  for  the 
experimentally  observed  minimum  chamber  length  for  instability  (  for 
chamber  lengths  less  than  the  minimum  the  combustion  is  stable)  Zucrow 
and  Osborn  postulated  that  the  net  energy  supplied  to  the  disturbance  de¬ 
pends  on  the  period  of  the  oscillation  and  hence  on  the  chamber  length. 
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For  short  chambers  the  period  of  the  oscillation  is  small,  hence  only  a 
small  amount  of  propellant  will  be  injected  during  each  cycle.  Thus  as 
the  chamber  length  decreases  the  energy  released  per  cycle  decreases,  and 
for  some  minimum  chamber  length  the  energy  addition  will  not  be  sufficient 
to  amplify  a  disturbance  (Implicit  in  this  argument  is  the  assumption  that 
the  energy  lost  per  cycle,  due  to  dissipation,  flow  work  and  mean  nozzle 
flow,  is  essentially  independent  of  the  period  of  the  oscillation). 

The  instability  mechanism  suggested  by  Zucrow  and  Osborn  was 
reviewed  by  Crocco  (32).  Crocco  pointed  out  that  the  energy  lost  per 
cycle  is  not  independent  of  the  period  of  the  oscillation,  but,  in  fact, 
decreases  with  decreasing  period.  Thus  the  balance  between  the  energy 
supplied  and  the  energy  lost  per  cycle  is  relatively  insensitive  to  the 
period  of  the  oscillation.  In  addition  Crocco  pointed  out  that  if  the 
Zucrow-Osborn  argument  were  carried  to  its  logical  conclusion  there  would 
be  no  upper  length  limit  for  stability-  ie.  no  maximum  length  above  which 
the  combustion  would  be  stable.  Crocco  noted  that  although  this  con¬ 
clusion  seemed  to  agree  with  available  gas  rocket  data  it  was  inconsis¬ 
tent  with  accepted  "time-lag"  theories  of  combustion  instability.  If,  as 
Zucrow  and  Osborn  postulated,  the  driving  mechanism  is  a  chemical  kinetic 
one  then  the  relevant  time-lag  should  be  the  chemical  reaction  time-  ie. 
the  time  required  to  convert  gaseous  reactants  into  gaseous  products.  For 


combustion  instability  to  occur  there  must  be  a  coupling  of  the  pressure 
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os  dilations  and  the  combustion  process.  For  this  coupling  to  occur  the 
ratio  of  the  chemical  reaction  time  and  the  period  of  the  oscillation 
must  lie  in  certain  ranges.  Since  the  period  of  the  oscillation  depends 
on  the  chamber  length  it  follows  that  there  must  be  a  lower  and  an  upper 
critical  length  for  each  mode  of  longitudinal  oscillation.  Below  or 
above  these  critical  lengths  combustion  will  be  stable  with  respect  to 
a  particular  mode  (Although  it  may  still  be  unstable  foi  -ther  modes). 

In  view  of  the  above  results  Crocco  concluded  that  the  instability 
mechanism  suggested  by  Zucrow  and  Osborn  was  incorrect. 

Since  it  appeared  that  the  driving  mechanism  for  gas-phase 
instability  was  not  related  to  chemical  kinetics,  a  different  driving 
mechanism,  based  on  a  coupling  between  oscillations  in  the  chamber  and 
the  heat  transfer  to  the  injector,  was  suggested  by  Bortzmeyer  and 
Crocco  (26).  In  this  driving  mechanism  it  is  assumed  that  the  gaseous 
propellants  are  in  thermal  equilibrium  with  the  injector-  ie.  the 
injection  temperature  of  the  gases  is  equal  to  the  temperature  of  the 
face  of  the  injector.  If  a  small  disturbance  occurs  near  the  injector 
face  the  heat  transfer  to  the  injector  will  increase  (due  to  the  increased 
gas  temperature  associated  with  the  disturbance).  At  some  later  time 
(determined  by  the  "inertia"  of  the  transfer  processes)  the  temperature 
of  the  injector  face  (and  thus  the  injection  temperature  of  the  gases) 
will  increase.  Due  to  the  increased  injection  temperature  (and  after  a 
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proper  time-lag)  the  temperature  of  the  burned  gases  increases.  Note 
that  as  the  temperature  of  he  injector  face  increases  the  heat  transfer 
to  the  injector  decreases  (since  the  difference  between  gas  temperature 
and  injector  temperature  decreases).  Bortzmeyer  and  Crocco  argued  that 
if  the  time-lags  associated  with  the  transfer  processes  were  properly 
related  to  the  period  of  the  oscillation,  unstable  combustion  would  occur. 
Using  this  drivirg  mechanism  Bortzmeyer  and  Crocco  carried  out  a  theo¬ 
retical  analysis  to  obtain  the  stability  limits  for  the  longitudinal 
mode  of  gas-phase  combustion  instability.  The  results  of  this  analysis 
agreed  in  some  respects  with  the  experimental  data  obtained  by  Zucrow 
and  Osborn  (31)  and  Pelmas  (27).  The  analysis  predicted  a  minimum 
chamber  length  for  instability  and  no  upper  length  limit  for  stability. 

It  also  predicted  that  Ci  u3tion  in  the  gas  rocket  tends  to  be  stable 
around  the  stoichiometric  mixture  ratio.  This  prediction  agreed  with 
the  results  of  Pelmas  but  not  with  the  results  of  Zucrow  and  Osborn 
(who  found  unstable  combustion  around  the  stoichiometric  mixture  ratio). 
The  stability  limits  obtained  from  the  analysis  are  strongly  dependent 
on  two  heat  transfer  parameters-  the  heat  transfer  coefficient  of  the 
combustion  gases  and  the  thermal  conductivity  of  the  injector.  An  ex¬ 
periment  in  which  the  heat  transfer  coefficient  of  the  combustion  gases 
(for  a  given  fuel  and  oxidizer)  was  altered  was  carried  out  by  Schob  (29). 
This  experiment  showed  effectively  no  dependence  of  the  stability  limits 
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on  the  heat  transfer  coefficient.  Experiments  using  geometrically 
similar  injectors  vith  and  vlthout  insulating  coatings  (28)  showed 
no  dependence  of  the  stability  limits  on  the  thermal  conductivity  of 
the  injector.  Thus  experiments  indicated  that  beat  transfer  to  the 
injector  has  little  or  no  effect  on  gas-phase  combustion  instability. 

Sirignano  and  Croce o  (23)  then  carried  out  a  theoretical 
analysis  of  the  longitudinal  mode  of  high-frequency  combustion  insta¬ 
bility  in  which  the  driving  mechanism  is  based  on  chemical  kinetics. 

The  mechanism  is  omevhat  similar  to  the  one  used  by  Zucrow  and  Osborn 
(31)  to  qualitatively  explain  their  experimental  results  in  that  in¬ 
stability  is  produced  by  a  coupling  of  the  pressure  oscillation  with  the 
chemical  reaction.  The  theory  assumes  that  the  response  time  of  the 
combustion  process  is  negligible  compared  to  the  period  of  the  oscillation. 
Hence,  no  phase  or  time-lag  exists  between  energy  addition  from  combus¬ 
tion  and  pressure.  The  analysis  provides  information  both  on  the  stability 
limits  and  on  the  instability  waveform.  The  theory  is  dicussed  in  some 
detail  in  the  following  sections  since  it  seems  particularly  relevant  to 
gas-phase  combustion  instability  and  forms  the  basis  for  much  of  the  ex¬ 
perimental  work  discussed  in  Chapter  4. 

The  explanation  for  some  of  the  experimentally  observed  pro¬ 
perties  of  the  longitudinal  mode  of  gas-phase  Instability  (in  particular 
the  existence  of  harmonic  modes)  cannot  be  found  in  the  Sirignano-Crocco 


model.  In  an  effort  to  account  for  the  presence  of  harmonics  Crocco  has 


recently  advanced  a  theory,  based  on  a  chemical  kinetics  driving  mech¬ 
anism,  which  assumes  a  combustion  zone  of  finite  thickness  (33).  Ana¬ 
lytical  and  experimental  studies  based  on  this  new  model  are  currently 
in  progress  at  Princeton. 

Culick  has  carried  out  a  theoretical  investigation  of  high- 
frequency  combustion  instability  for  the  general  case  of  a  three- 
dimensional  flow  field  with  distributed  combustion  (21).  A  solution 
of  the  linearized  equations  gives  an  expression  fo**  the  stability  limits 
in  which  the  effects  of  combustion,  nozzle  and  mean  flow  are  separated. 
Under  suitable  simplifications  it  should  be  possible  to  reduce  the 
equation  for  the  stability  limits  to  a  form  which  is  applicable  to  the 
longitudinal  mode  of  gas-phase  instability.  Since  the  Culick  analysis 
provides  less  information  about  the  instability  phenomenon  (stability 
limits  only)  than  the  Sirignano-Crocco  analysis,  no  attempt  was  made  in 
the  present  study  to  correlate  Culick' s  results  with  experimental 
observations. 

Some  comments  on  the  results  of  the  various  studies  discussed 
in  this  section  are  given  in  Chapter  5* 

B.  The  Sirignano-Crocco  Instability  Model 


Sirignano  and  Crocco  (23)  have  carried  out  a  theoretical  in- 
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vestigation  of  the  longitudinal  mode  of  high-frequency  combustion  in¬ 
stability  based  on  a  chemical  kinetics  driving  mechanism.  The  driving 
mechanism  assumes  that  instability  is  produced  by  a  coupling  of  the 
pressure  oscillation  with  the  chemical  reaction  and  that  there  is  no 
time-lag  between  energy  addition  from  combustion  and  pressure.  The 
model  on  which  the  analysis  is  based  closely  approximates  the  conditions 
in  a  rocket  burning  premixed  gaseous  propellants.  In  the  following 
sections  a  brief  discussion  of  the  theory  is  presented  to  serve  as  a 
background  for  the  experimental  work  discussed  in  Chapter  4.  Details 
of  the  theoretical  development  can  be  found  in  Reference  23  and 
Appendix  C. 

The  Sirignano-Crocco  analysis  differs  from  the  linearized 
approaches  to  the  instability  problem  (8,  17-19>  21)  in  that  it  seeks  a 
solution  to  the  non-linear  equations  of  motion-  ie.  a  solution  which  is 
valid  inside  a  region  of  unstable  combustion  and  not  Just  at  the  stabil¬ 
ity  limits.  The  principal  assumptions  on  which  the  analysis  is  based  are: 

(1)  The  response  time  of  the  combustion  process  is  negligible 
compared  to  the  period  of  the  oscillation. 

(2)  A  plane  shock  wave  propagates  longitudinally  in  the 
combustion  chamber  with  constant  period,  and  reflects  alternately  from 
the  nozzle  and  injector. 

(3)  The  combustion  chamber  is  sufficiently  long  so  that  the 
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configuration  may  be  approximated  by  the  limiting  case  of  a  zero-length 
nozzle  and  a  zero- thickness  combustion  zone. 

(4)  There  are  no  entropy  waves  in  the  flow  field. 

(5)  The  flow  is  one -dimensional  and  adiabatic. 

(6)  The  combustion  gases  are  calorically  perfect. 

(7)  The  cross-sectional  area  of  the  combustion  chamber  is 

constant. 


Because  of  the  presence  of  shock  waves  in  the  flow  field  it  is 
convenient  to  use  a  characteristic  co-ordinate  perturbation  technique  to 
obtain  solutions.  To  this  end  the  unsteady,  one-dimensional  equations  of 
motion  of  the  gas  in  the  chamber  are  transformed  into  characteristic  co¬ 
ordinates,  and  the  dependent  variables  (which  include  x,t  in  the  char¬ 
acteristic  co-ordinate  system)  are  expanded  in  terms  of  an  as  yet  un¬ 
defined  amplitude  parameter  £  .  The  boundary  condition  at  the  combustion 
zone  is  defined  in  terms  of  an  energy  release  parameter  r.  The  boundary 
condition  at  the  nozzle  is  defined  in  terms  of  a  nozzle  admittance 
parameter.  It  is  noted  that  the  Sirignano-Crocco  analysis  differs  from 
the  standard  approaches  to  the  non-linear  instability  problem  (20 )  in 
that  it  seeks  a  "stable"  periodic  solution  rather  than  an  indication  of 
the  stability  of  a  given  steady-state  operation  to  known  perturbations. 
Hence,  the  initial  condition  is  replaced  by  a  cyclic  condition  (which  in 
effect  states  that  the  oscillation  is  periodic).  "Stable"  periodic 
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solutions  to  the  equations,  compatible  with  the  boundary  conditions  and 
the  cyclic  condition,  are  obtained  only  after  considerable  labor  (Phe¬ 
nomenologically  a  "stable"  periodic  solution  exists  when  the  energy 
supplied  to  the  wave  by  the  combustion  process  Just  balances  the  energy 
dissipated  by  the  wave  due  to  its  reflection  from  the  nozzle  and  its 
motion).  The  solutions  give  the  strength  of  the  shock  wave  (in  terms 
of  the  changes  in  pressure  and  gas  velocity  across  the  wave),  the  wave 
velocity,  the  period  of  the  oscillation  and  the  time  decay  (for  fixed 
location)  in  pressure  and  gas  velocity  behind  the  wave  as  functions  of 
three  important  instability  parameters-  f  ,  r  (Appendix  C).  If 

the  functional  form  of  the  combustion  law  (r)  is  known  then  it  is 
possible  to  obtain  ^  and  \  .  Knowing  ^  \  and  r  the  waveform  of 

the  instability  can  be  determined.  Alternatively,  if  the  waveform  of 
the  instability  is  known  (through  experiments)  then  it  is  possible,  in 
theory,  to  work  backwards  to  obtain  the  appropriate  functional  form  of 
the  combustion  law.  In  the  following  sections  some  theoretical  results 
obtained  for  a  choice  of  r  particularly  appropriate  for  gas  rockets  are 
presented. 


C.  Theoretical  Results 

For  any  reasonable  r  the  waveform  obtained  from  the  theory 
consists  of  shock  discontinuities  followed  by  exponential  time  decay 
in  pressure  and  velocity.  The  shock  wave  oscillates  with  approximately 
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the  fundamental  acoustic  frequency  of  the  combustion  cljamber. 

In  order  to  determine  some  explicit  information  about  the 
instability  waveform  Sirignano  and  Crocco  assumed  a  simplified  combus¬ 
tion  law  of  the  Arrhenius  type 

r  =  Be»p{-£}  (3-1) 

where  E=  an  activation  energy,  T=  combustion  temperature,  B-  pre-expo¬ 
nential  factor  (independent  of  T)  and  R=  universal  gas  consuint.  Ad¬ 
mittedly  this  expression  may  be  too  simple  to  be  realistic,  but  it  does 
contain  some  of  the  essential  features  of  gas-phase  kinetics  and  provides 
theoretical  results  which  are  readily  compared  with  the  experimental  data. 
With  this  choice  for  r  the  functional  forms  of  £  and  ^  are  (23) 

€  (3-2) 


(3-3) 

The  parameter  £  can  be  considered  an  amplitude  parameter,  roughly  pro- 
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portional  to  the  strength  of  the  instability  (Appendix  C).  The  parameter 
\  can  be  considered  a  shape  parameter  giving  the  nature  of  the  exponen¬ 
tial  decay  in  pressure  and  velocity  behind  the  wave  (Appendix  C). 


D.  Discussion  of  the  Assumptions 

In  this  section  some  of  the  assumptions  on  which  the  theory  is 
based  are  reviewed,  with  particular  emphasis  being  placed  on  a  comparison 
of  the  assumptions  with  conditions  in  the  gas  rocket. 

Assumption  1.  The  response  time  of  the  combustion  process  is 
negligible  compared  to  the  period  of  the  oscillation. 

The  response  time  of  the  combustion  process  is  the  time  re¬ 
quired  for  the  reaction  rate  to  adjust  to  a  change  in  conditions  (pressure, 
temperature )- 

I  I 
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where  k=  specific  rate  constant,  [reactant  ij  =  concentration  of  reactant 
i  and  n^=  order  of  the  reaction  with  respect  to  reactant  i.  For  the  con¬ 
ditions  existing  in  the  gas  rocket  typical  values  for  t  are  of  the 

response 

order  of  10  ^  seconds.  Topical  values  for  the  period  of  the  oscillation 

-2  -3 

are  10  to  10  seconds.  Thus  for  rocket  motors  burning  premixed  gaseous 
propellants  the  assumption  of  no  time-lag  between  energy  addition  by 


I  • 


combustion  and  the  oscillation  is  a  valid  one  (This  is  not  true  for 
liquid  propellant  motors  where  energy  addition  depends  on  factors  other 
than  gas-phase  kinetics). 

Assumption  2.  A  plane  shock  wave  propagates  longitudinally  in 
the  combustion  chamber. 

Experimental  studies  of  the  longitudinal  mode  of  high-frequency 
combustion  instability  in  liquid  propellant  rockets  (9-1 2,  14)  and  in  gas 
rockets  (27,  29-31,  3*0  have  shown  that  the  instability  frequently  appears 
as  a  shock  wave  propagating  longitudinally  Jn  the  combustion  chamber. 

This  observation  is  confirmed  by  the  present  experimental  study  (Chapter  4). 

Assumption  3.  The  nozzle  length  and  combustion  zone  thickness 
are  negligible  compared  to  the  combustion  chamber  length. 

The  combustion  zone  thickness  and  nozzle  length  in  the  gas 
rocket  are  generally  very  much  less  than  the  length  of  the  combustion 
chamber  (See  Chapter  4).  Hence,  Assumption  3  (which  is  necessary  to 
make  the  problem  mathematically  tenable)  is  valid,  to  a  first  approxi¬ 
mation,  for  the  gas  rocket. 

In  general  a  combustion  zone  and  nozzle  of  finite  size  have  a 
stabilizing  effect  on  the  combustion  process.  The  stabilizing  effect  of 
a  finite  combustion  zone  results  from  the  non-optimum  addition  of  energy 
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to  the  combustion  oscillation-  ie.  energy  addition  over  a  finite  length 
rather  than  at  a  pressure  antinode.  It  is  conceivable  that  for  suffi¬ 
ciently  short  chambers  (ie.  for  pressure  oscillations  with  small  wave 
length)  the  combustion  oscillations  would  not  be  amplified.  In  fact  this 
argument  is  frequently  offered  as  an  explanation  for  the  experimentally 
observed  minimum  chamber  length  for  longitudinal  instability  (32).  As 
the  length  of  the  combustion  chamber  increases  the  energy  addition 
approaches  the  optimum  and  the  stabilizing  effect  decreases. 

Any  stabilizing  effect  of  the  nozzle  results  from  variations 
in  the  amplitude  and  phase  of  pressure  and  velocity  produced  by  the 
interaction  of  the  oscillation  with  the  nozzle.  In  general  these  ampli¬ 
tude  and  phase  variations  decrease  with  decreasing  frequency  (longer 
chambers);  hence,  as  the  chamber  length  increases  the  stabilizing  effect 
of  the  nozzle  decreases.  An  experimental  study  of  the  effects  of  the 
combustion  zone  and  nozzle  on  the  instability  phenomenon  is  discussed 
in  Chapter  4. 

Assumption  4.  No  entropy  waves  exist  in  the  flow  field. 

In  a  rocket  motor  burning  premixed  gaseous  propellants  entropy 
waves  can  be  produced  by  localized  pressure  variations  (such  as  accompany 
unstable  combustion)  in  the  combustion  zone.  Adjacent  propellant  elements 
will  have  different  entropies;  the  mean  downstream  motion  of  these 
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elements  is  termed  an  entropy  wave.  An  entropy  wave  can  have  two  possible 
effects  on  the  instability  phenomenon.  The  first  effect  is  the  appear¬ 
ance  of  intermediate-frequency  instability  (discussed  in  Chapter  l).  The 
second  effect  results  from  the  interaction  of  the  entropy  wave  with  the 
pressure  waves  associated  with  high-frequency  instability.  Due  to  the 
highly  turbulent  flow  in  the  combustion  chamber  of  the  gas  rocket  it  is 
reasonable  to  assume  that  entropy  waves  are  "erased"  due  to  turbulent 
mixing.  If  this  is  true  then  the  effects  of  the  entropy  wave  can  be 
neglected.  In  the  experiments  of  Chapter  k  no  effects  attributable  to 
entropy  waves  were  observed. 

Assumption  5.  The  flow  is  one -dimensional  and  adiabatic. 

Due  to  the  high  turbulence  level  downstream  from  the  combustion 
zone  it  is  reasonable  to  assume  that  the  mean  flow  properties  are  con¬ 
stant  over  a  given  cross-section.  Tbe  nature  of  the  flow  immediately 
downstream  from  the  injector  and  in  the  combustion  zone  depends  on  the 
injection  pattern.  The  effects  of  non-uniform  flow  in  the  combustion 
zone  on  the  instability  phenomenon  are  not  known  a  priori.  Hence,  an 
experimental  study  of  the  efft  ts  of  the  injection  pattern  on  the  in¬ 
stability  has  been  made  (Chapter  k) . 

The  flow  in  the  gas  rocket  is  definitely  non- adiabatic.  Heat 


transfer  to  the  walls  of  the  combustion  chamber  is  relatively  large  for 


-29- 


stable  combustion,  and  when  the  combustion  is  unstable  the  heat  transfer 
increases  significantly.  Heat  transfer  has  a  stabilizing  effect  on  the 
combustion  process  (as  do  all  dissipative  phenomena). 

Assumptions  in  Equation  (3-1) • 

Two  important  assumptions  are  embodied  in  the  energy  release 
law  (3-1)*  First,  the  energy  release  rate  is  assumed  to  be  independent 
of  pressure  or  species  concentration.  A  more  rigorous  expression  for 
the  energy  release  rate  would  be 

A 

r  *qr 

where  q  is  a  heat  of  reaction  and  $  is  the  chemical  reaction  rate  and  is 
written 

r  *  b(T)  n  [reactant  i]  '  exp  {- 

where  b(T)=  a  temperature  dependent  pre- exponential  factor.  Expressing 
the  concentration  of  the  ith  reactant  in  terras  of  the  mole  fraction  (x^ 

and  the  total  pressure  the  energy  release  law  becomes 

Ui 

r=qb(T)[|f]  nXi"'  exp  {-It} 


-30- 


In  employing  an  ene.gy  release  law  of  this  type  Sirignano  and  Crocco 

found  that  for  reasonable  values  of  £  n..  (the  overall  order  of  the 

l 

reaction)  the  combustion  was  predicted  always  to  be  unstable.  This 
conclusion  is  contrary  to  experimental  observations  (Chapter  4).  That 
the  more  rigorous  rate  expression  predicts  a  result  not  in  agreement 
with  experimental  data  is  due  primarily  to  the  fact  that  the  isentropic 
assumption  used  to  relate  P  and  T  in  the  analysis  of  the  combustion  zone 
boundary  condition  is  not  valid.  The  constant  entropy  assumption  is 
violated  due  to  the  energy  release  and  diffusion  in  the  combustion  zone. 
Hence  to  obtain  non-trivial  theoretical  results  it  is  necessary  to  assume 
that  the  combustion  process  in  the  gas  rocket  is  approximated  sufficiently 
well  by  equation  (3-1 ).  As  will  be  seen  in  Chapter  k  simplified  rate 
expressions,  6uch  as  (3-l),  provide  theoretical  results  which  agree  in 
many  respects  with  the  experimental  observations. 

A  direct  result  of  the  form  of  the  energy  release  law  (3-l)  is 
that  the  mean  energy  release  rate  (“)  can  increase  without  limitation  as 
the  mean  combustion  temperature  (T)  increases.  In  the  gas  rocket  the 
mean  energy  release  rate  is  the  product  of  the  propellant  injection  rate 
(m)  and  a  heat  of  reaction  (q), 


P  =  m  q 
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During  a  given  run  m  and  the  mixture  ratio  (and  hence  q)  are  constant. 

It  follows  that  the  mean  energy  release  rate  is  constant  (and  essentially 
independent  of  the  mean  combustion  temperature } .  If  the  mean  combustion 
temperature  were  increased  (eg.  by  heat  audition  through  the  walls  of 
the  combustion  chamber)  equation  ( 3-1 )  would  predict  a  corresponding 
increase  in  the  mean  energy  release  rate.  In  the  gas  rocket  the  in¬ 
creased  combustion  temperature  would  produce  a  temporary  increase  in  the 
Instantaneous  energy  release  rate-  ie.  an  increase  in  the  propellant 
consumption  rate.  Since  the  propellant  consumption  rate  is  now  greater 
than  the  propellant  injection  rate  the  increased  energy  release  rate 
cannot  persist.  After  some  time  the  energy  release  rate  must  decrease 
so  that  in  the  mean  the  propellant  consumption  rate  is  equal  to  the 
propellant  injection  rate.  Thut  it  is  seen  that  equation  ( 3-1 )  neglects 
the  fact  that  the  gas  rocket  is  an  energy-limited  system-  ie.  that  the 
energy  release  rate  is  limited  by  the  propellant  injection  rate. 

A  further  discussion  of  the  effects  of  deviations  from  the 
assumptions  on  the  instability  phenomenon  is  presented  in  Chapter  4. 

E.  Discussion  of  Theoretical  Results  and  Some  Critical 
Experiments 

A  careful  study  of  the  results  of  the  Sirignano-Crocco  analysis 
(Section  C  and  Appendix  C)  suggests  several  experiments  which  can  be 


used  to  test  the  validity  of  the  proposed  instability  driving  mechanism 


and  combustion  lav. 

The  theory  predicts  that  the  instability  exists  as  a  plane 
shock  vave  oscillating  longitudinally  in  the  combustion  chamber  with 
approximately  the  fundamental  acoustic  frequency  of  the  chamber. 
According  to  the  theory  the  time-decay  in  pressure  behind  the  vave  (at 
fixed  location)  is  exponential  and  depends  on  the  functional  form  of  the 
combustion  lav.  To  obtain  data  on  the  instability  waveform  an  experi¬ 
ment  can  be  performed  in  which  pressure-tine  histograms  are  obtained 
during  unstable  combustion  at  various  locations  in  the  chamber  (Chapter 

h). 


From  equations  (3-2),  (3-3) >  (C-3)  and  (C-U)  it  is  seen  that 
the  equation  for  the  stability  limits  (a  stability  limit  is  defined  by 
A  p/p  -  0)  is  obtained  by  setting  £  *  0  in  equation  (3-2). 

(L\  s  (3y- 1)  ,  , 
\rT/..  2  (r~  i)  3' 

limn 

Thus  the  stability  limits  are  defined  in  terms  of  a  critical  value  for 
tne  rate  parameter  E/rT.  For  a  given  propellant  combination  (and  over 
a  moderate  range  of  mixture  ratios)  it  is  reasonable  to  assume  that  E 


is  constant.  Iftus  the  stability  limits  are  defined  in  terms  of  a  critical 
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combustion  temperature,  From  e  pat  ion  (3-2)  it  is  seen  that 

whenever  T  <  Tcrl^  the  amplitude  parameter  «  >0;  hence,  from  equations 

(C-3)  and  (C-4),  £  P/P  >  0-  ie.  the  combustion  is  unstable.  Whenever 

T  >  Tcrit,  «  <  0  and  A  <  0-  ie.  the  combustion  is  stable.  The 

important  conclusion  to  be  drawn  from  this  dicussion  is  that  the  tendency 

for  instability  is  greater  at  lower  combustion  temperatures.  It  is  noted 

that  the  combustion  temperature  depends  on  the  equivalence  ratio  (  $  ), 

is  a  maximum  at  stoichiometric  (  <J>  «  l)  and  decreases  as  the  values 

for  *  vary  from  $  -  1.  If  the  maximum  combustion  temperature  for 

the  propellant  combination  is  greeter  than  T  ,,  then  two  values  for 

crix 

*  *  ^.<25 

♦lUnlt  'xl6t-  *  limit  «  1  and  •  limit  »  U  Thl6  10  1"t'n>r«ted 
tc  mean  that  in  the  $  -L  or  <J>  -T  plane  two  regions  of  unstable  combus¬ 
tion  will  be  found,  one  on  either  side  of  stoichiometric.  According  to 
the  theory  once  a  stability  limit  has  been  crossed  (proceeding  away  from 

♦  «  l)  unstable  combustion  persists  until  the  combustion  limit  is 
reached.  When  the  maximum  combustion  temperature  is  less  than 
unstable  combustion  will  occur  over  the  entire  range  of  combustible 
equivalence  ratios.  From  the  above  considerations  it  is  seen  that  the 
theory  provides  a  criterion  for  the  number  and  location  of  unstable 
regions . 


From  equation  (3-M  it  is  Been  that  for  propellant  combinations 
characterized  by  two  unstable  regions  (T^^  >  Tcrit)  it  is  possible  to 


collapse  the  two  regions  into  a  single  region  if  the  overall  range  of 
combustion  temperatures  can  be  reduced  without  altering  E  (eg.  by  the 
addition  of  an  inert  diluent).  The  two  unstable  regions  should  Just 
coalesce  into  a  single  region  when  TWWY«  ^crit* 


To  test  the  concept  of  a  critical  value  for  E/RT  at  the 

stability  limits  two  experiments  can  be  carried  out.  In  the  first 

experiment  the  regions  of  unstable  combustion  for  several  propellant 

combinations  are  located  in  the  $  -L  and  $  -P  planes.  For  those 

combinations  which  exhibit  two  unstable  regions  a  comparison  is  made 

between  the  experimental  values  for  (E/RT)limlt  and  the  values  given 

by  equation  (3-M-  For  those  propellant  combinations  which  exhibit 

one  unstable  region  the  theoretical  criterion  (E/RT), ,  ,  <  (e/rT) 

limit  stoic. 

is  tested. 


The  second  experiment  involves  one  of  the  prooellant  combina¬ 
tions  which  exhibits  two  unstable  regions.  In  this  experiment  the 
overall  range  of  combustion  temperatures  is  reduced  (without  altering 
E)  through  the  addition  of  an  inert  diluent,  and  the  stability  limits 
for  various  diluent  concentrations  obtained.  According  to  equation 
(3-M»  (e/FT) ^ should  remain  essentially  constant  (a  slight  varia¬ 
tion  in  Y  is  produced  by  the  addition  of  a  diluent).  Since  the  combus¬ 


tion  temperatures  are  lower  the  theory  indicates  a  shift  in  the  stability 
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limits  toward  $  =  1  (and  relatively  higher  combustion  temperatures) 
as  the  amount  of  diluent  is  increased. 

From  equation  (3-M  it  is  seen  that  the  stability  limits 
depend  on  an  overall  activation  energy  characteristic  of  the  propellant 
combination.  To  determine  the  dependence  of  the  stability  limits  on  E 
an  experiment  involving  a  propellant  combination  which  exhibits  two 
unstable  regions  can  be  performed.  In  this  experiment  E,  but  not  T, 
is  altered  by  the  addition  of  trace  amounts  of  a  catalyzing  species. 

In  accordance  with  equation  (3-M»  a  shift  in  stability  limits  should 
be  observed  with  the  addition  of  the  catalyzing  species. 

Equations  (3-2),  (3-3) ,  ( C- 3)  and  (C-U)  give  the  dependence 
of  the  strength  of  the  instability  ( A  P/P)  on  the  rate  parameter  (E/RT) 
and  the  mean  Mach  number  of  the  combustion  gases  (M).  A  typical  plot 
of  the  theoretical  variation  of  A  P/P  with  E/FT  for  constant  M  is  shown 
in  Figure  l4.  It  is  noted  that  A  P/P  is  a  monotonically  increasing 
function  of  E/RT.  For  a  given  propellant  combination  E  is  approximately 
constant,  hence  A  P/P  is  a  monotonically  decreasing  function  of  T. 
Recalling  the  dependence  of  f  on  ♦  ,  it  is  seen  that  A  P/P  increases 

as  ♦  varies  away  from  stoichiometric.  Thus  for.  propellant  combina¬ 
tions  which  exhibit  two  unstable  regions  it  is  to  be  expected  that  the 
instability  strength  increases  as  an  equivalence  ratio  traverse  is  made 


across  an  unstable  region  (from  the  stability  limit  to  the  combustion 
limit);  the  strength  is  a  maximum  at  the  combustion  limits.  For 
propellant  combinations  which  exhibit  one  unstable  region  it  is  to  be 
expected  that  the  strength  is  a  minimum  at  stoichiometric  and  a  maximum 
at  the  combustion  limits.  To  determine  the  dependence  of  A  P/P  on 
E/RT  for  constant  M  an  experiment  can  be  carried  out  in  which  a  P/P 
is  measured  as  equivalence  ratio  traverses  are  made  across  regions  of 
unstable  combustion. 

A  typical  plot  of  the  theoretical  variation  of  A  P/P  with  M 
for  constant  E/RT  i6  shown  in  Figure  15.  It  is  seen  that  A  P/P  is  a 
monotonically  increasing  function  of  M.  It  is  noted  that  for  given 
propellants  M  depends  primarily  on  the  throat  area  of  the  exhaust 
nozzle  (for  choked  nozzle  flow).  To  determine  the  dependence  of  A  P/P 
on  M  for  constant  E/rT  an  experiment  can  be  performed  in  which  A  P /P 
is  measured  for  fixed  E/RT  (ie.  for  fixed  •  )  using  exhaust  nozzles 

with  different  throat  areas. 

As  noted  in  Section  D  the  instability  phenomenon  depends  to 
a  certain  extent  on  the  injection  pattern,  nozzle  geometry  and  combus¬ 
tion  chamber  length.  An  experimental  study  of  these  effects  has  been 
carried  out  and  is  discussed  in  Chapter  4. 


A  detailed  description  of  the  experiments  outlined  in  this 
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section  together  with  a  comparison  of  the  experimental  results  with 


the  theory  is  presented  in  Chapter  4. 


-38- 


CHAPTER  4 .  EXPERIMENTAL  CONSIDERATIONS 

A.  General  Comments  on  the  Operation  of  the  Gas  Rocket 

A  typical  gas  rocket  run  consists  of  three  distinct  phases- 
starting,  steady-state  and  data  recording.  During  the  starting  phase 
combustion  is  established  with  reduced  propellant  flow  rates.  Ihis 
procedure  serves  to  reduce  the  severity  of  flashback  (should  it  occur) 
and  to  prevent  a  "hard"  start.  The  details  of  the  starting  procedure 
depend  on  the  individual  experiment  and  will  be  discussed  in  the  sections 
describing  the  experiments.  Steady-state  operation  is  achieved  by 
gradually  increasing  the  propellant  flow  rates  to  t.heir  steady-state 
values  (determined  by  mixture  ratio  and  combustion  pressure  require¬ 
ments).  Steady-state  operation  is  characterized  by  constant  propellant 
flow  rates  and  constant  mean  combustion  pressure.  The  duration  of 
steady-state  operation  is  on  the  order  of  20  to  30  seconds  so  that 
transients  associated  with  starting  and  with  transition  to  steady-state 
operation  have  sufficient  time  to  decay.  Instability  data  are  recorded 
during  the  final  10  to  15  seconds  of  the  steady-state  period. 

B.  Waveform  of  Longitudinal  Combustion  Instability 

T^c  instability  waveform  obtained  from  the  Sirignano-Crocco 
analysis  consists  of  shock  discontinuities  followed  by  exponential 
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decay  in  time  (for  fixed  location)  of  pressure  and  velocity.  The  shock 
wave  oscillates  with  approximately  the  fundamental  longitudinal  acoustic 
frequency  of  the  combustion  chamber.  A  typical  theoretical  pressure¬ 
time  history  is  shown  in  Figure  l6a.  The  results  are  for  an  observer 
located  at  x/L  =  0.l6  (ie.  near  the  injector). 

In  the  present  study  pressure- time  histories  during  unstable 
combustion  were  recorded  at  several  stations  in  the  combustion  chamber 
of  the  gas  rocket.  Figure  l6b  is  a  typical  pressure-time  history 
during  unstable  combustion  at  a  station  Just  downstream  from  the  in¬ 
jector.  The  first  pressure  rise  corresponds  to  a  wave  traveling  up¬ 
stream  toward  the  injector.  The  second  pressure  rise  corresponds  to 
the  wave  after  reflection  from  the  injector  (traveling  downstream 
toward  the  nozzle).  A  comparison  of  the  theoretical  and  experimental 
pressure-time  histories  shows  that  they  are  qualitatively  the  same. 

From  pressure-time  histories  such  as  Figure  l6b  it  has  been  determined 
that  the  instability  occurs  as  a  plane  shock  wave  oscillating  with 
approximately  the  fundamental  longitudinal  acoustic  frequency  of  the 
combustion  chamber  with  exponential  time-decay  of  pressure  behind  the 
wave.  For  certain  conditions  of  operation  other  forms  of  combustion 
oscillations  have  been  observed.  These  are  discussed  in  Section  F. 
Waveform  of  the  type  shown  in  Figure  l6b  have  been  observed  in  other 
studies  of  combustion  instability  in  gas  rockets  (27,  29,  31,  3M  • 


C.  Dependence  of  Stability  Limits  on  E/FT 


In  Chapter  3  (Section  E)  it  was  noted  that  the  stability 
limits  are  defined  in  terms  of  a  critical  value  for  the  rate  parameter 
(E/RT).  Several  experiments  which  serve  to  test  this  conclusion  were 
discussed.  In  each  of  these  experiments  the  dependence  of  the  stability 
limits  on  one  of  the  two  variables  in  the  rate  parameter  (E  and  T)  is 
determined  holding  the  other  variable  constant. 

1.  Dependence  of  Stability  Limits  on  T  (E  =  constant) 

To  determine  the  dependence  of  the  stability  limits  on  T  the 
regions  of  unstable  combustion  for  two  propellant  combinations  with 
various  concentrations  of  an  inert  diluent  were  located  in  the  $  -L 
plane.  The  two  propellant  combinations  used  in  this  experiment  were 
hydrogen-air  and  methane- oxygen.  The  inert  diluent  was  nitrogen.  The 
addition  of  nitrogen  to  a  given  propellant  combination  reduces  the 
overall  range  of  combustion  temperatures  but  has  little  or  no  effect 
on  the  overall  activation  energy  (35)  •  Thus  the  temperature  dependence 
of  the  stability  limits  can  be  obtained  by  determining  the  dependence 
of  the  limits  on  d.luent  concentration.  Quantitative  information  on  T 
at  a  stability  Unit  was  derived  from  the  axial  distribution  of  mean 
gas  temperature  i;i  the  combustion  zone  at  the  stability  limit. 
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It  is  noted  that  the  diluent  vas  supplied  to  the  rocket  motor 
separately  from  the  fuel  and  oxidizer  (in  the  flow  schematic  of  Figure 
1  the  diluent  is  the  X  gas).  During  the  starting  phase  combustion  was 
established  using  the  fuel  and  oxidizer  alone  (no  diluent).  Steady- 
state  operation  was  then  established  by  increasing  the  flow  rates  of 
the  fuel,  oxidizer  and  diluent  to  their  steady-state  values.  Ibis 
technique  is  somewhat  different  from  the  one  used  in  previous  inves¬ 
tigations  of  gas-phase  instability  (27-29).  In  these  studies  the 
desired  diluent  concentrations  were  obtained  using  gas  cylinders  in 
which  the  diluent  and  oxidizer  were  premixed  to  specification  by  the 
supplier.  The  principal  advantage  of  adding  the  diluent  separately 
is  that  instability  data  can  be  obtained  for  diluent  concentrations 
which  could  not  have  been  ignited  (using  a  spark  ignition  system) 
had  the  diluent  and  oxidizer  been  premixed.  In  addition  the  use  of 
separate  diluent  and  oxidizer  systems  increases  the  flexibility  of  the 
propellant  feed  system  in  that  variations  in  diluent  concentration  are 
easily  obtained.  In  using  this  technique  it  was  found  that  diluent 
concentrations  were  reproducible  to  within  t3 %  of  a  nominal  value. 

The  regions  of  unstable  combustion  in  the  $  -L  plan**  for 
hydrogen-air  and  hydrogen-dilute  air  are  shown  in  Figure  17.  In  this 


figure  the  diluent  concentration  is  defined  by 
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diluent  mass  flow  rate 

$  dilution  »  .. — . . .  . —  x  100. 

diluent  mass  flow  rate  +  air  mass  flow  rate 

The  unstable  regions  in  the  $  -L  plane  for  methane-dilute  oxygen  are 
shown  in  Figure  19.  The  data  of  Figures  1?  and  19  were  obtained  for 
P  *  7*8  *  0.3  atm  using  the  31-hole  showerhead  injector  and  the  plug 
nozzle.  Unstable  combustion  (characterized  by  periodic  oscillations 
of  the  combustion  pressure)  occurs  inside  the  shaded  regions. 

From  Figure  17  it  is  seen  that  for  a  given  diluent  concentra¬ 
tion  two  regions  of  unstable  combustion  exist,  one  on  either  side  of 
the  stoichiometric  mixture  ratio.  As  the  diluent  concentration  is  in¬ 
creased  the  unstable  regions  shift  toward  stoichiometric.  The  upper 
stability  limits  of  the  oxygen-rich  (  ♦  >  l)  unstable  regions  are  not 
clearly  defined-  ie.  there  is  no  sharp  transition  from  unstable  combus¬ 
tion  to  stable  '■’ombustion.  For  clarity  only  the  upper  limit  for  f^-air 
is  shown.  The  overall  extent  of  the  oxygen-rich  \nstable  regions  is 
only  slightly  dependent  on  the  diluent  concentration  since  as  the 
diluent  concentration  increases  both  the  upper  and  lower  stability  limits 
shift  toward  stoichiometric.  ITie  lower  stability  limits  for  the  fuel- 
rich  (  $<  l)  unstable  regions  are  in  reality  the  fuel-rich  "ombustion 
limits.  For  clarity  only  the  lower  limit  for  H^-air  is  shown.  The 
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overall  extent  of  the  fuel-rich  unstable  regions  is  only  slightly 
dependent  on  the  diluent  concentration.  The  results  of  Figure  17  for 
I^-air  are  in  general  agreement  with  the  observations  of  Pelmas  (27), 
Bertrand  (28)  and  Schob  (29).  Bertrand  (28)  has  observed  a  shift  in 
the  stability  limits  with  increasing  diluent  concentration. 

From  Figure  19  it  is  seen  that  for  (^-(O.UOp+O^N^)  two  regions 
of  unstable  combustion  exist,  one  on  either  side  of  stoichiometric. 

The  upper  limit  of  the  oxygen-rich  unstable  region  and  the  lower  limit 
of  the  fuel-rich  unstable  region  are  combustion  limits.  As  the  diluent 
concentration  i6  increased  the  unstable  regions  shift  toward  stoichio¬ 
metric.  For  CH^-air  the  fuel-rich  and  oxidizer-rich  unstable  regions 
have  coalesced  to  form  a  single  unstable  region  located  around 
stoichiometric.  The  upper  and  lower  boundaries  of  this  region  arc 
combustion  limits.  The  results  of  Figure  19  for  H^-air  are  in  general 
agreement  with  the  observations  of  Zucrow  and  Osborn  (24).  The  results 
for  CHj+-(0.40^+0.6N^))  are  in  general  agreement  with  the  observations  of 
Pelmas  (27) • 

A  quantitative  measure  of  was  obtained  from  the  axial 

distribution  of  mean  gas  temperature  in  the  combustion  zone  at  the 
stability  limits  (See  Figure  43).  The  details  of  the  measurement  of 
gas  temperatures  are  presented  In  Section  K.  Hie  gas  temperature  ^and 


thus  the  reaction  rate)  varies  throughout  the  combustion  zone.  For  the 
purpose  of  correlating  the  experimental  data  with  the  theoretical  re¬ 


sults  an  average  value  for  was  defined  by 

ft 

•‘■>{-4}  s  T  / #xp  {"  Rfo)  dx 


(4-1) 


In  this  definition  7  Is  that  temperature  which  would  give  an  average 
energy  release  rate  equal  to  the  local  energy  release  rate  integrated 
over  the  combustion  zone.  Two  experimental  inputs  are  required  in  (4-1)- 
the  thickness  of  the  combustion  zone  (  §  )  and  the  temperature  profile 
in  the  combustion  zone,  T(x).  The  thickness  of  the  combustion  zone  is 
taken  as  the  distance  from  the  injector  face  to  the  peak  combustion 
temperature.  Since  temperature  data  could  not  be  obtained  near  the  in¬ 
jector  (See  Section  K)  an  accurate  determination  of  T(x)  could  not  be 
made.  For  the  purpose  of  evaluating  Tqv  a  linear  temperature  profile 
was  assumed- 


8  ^inj  +  ^p#ok~^inj)  g 

where  T^nj  =  mean  injection  temperature  and  T^^  =  mean  peak  combustion 


temperature.  Using  this  profile  the  average  combustion  temperatures 


at  the  stability  limits  for  Hg-air  and  Hg-dilute  air  were  calculated  for 
several  choices  of  It  was  found  that  Tfty  was  approximately  equal 

to  Tpggfc  nearly  Independent  of  the  choice  of  .  The  results  of 
this  calculation  together  with  the  stoichiometric  adiabatic  combustion 
temperatures  (T^)  for  the  various  diluent  concentrations  are  tabulated 
in  Figure  18.  For  a  given  chamber  length  the  measured  combustion 
temperatures  at  the  stability  llmitr  nearest  stoichiometric  (^ay)  are 
approximately  the  same.  For  L  »  104  cm  -  llv05i25°K  for  all 
diluent  concentrations.  For  L  ■  53  cm  ^av  ■  1^35±^0C^‘  scatter 

in  the  temperature  data  can  be  explained  in  terms  of  the  experimental 
uncertainties  in  the  equivalence  ratio,  diluent  concentration  and 
thermocouple  measurement  (See  Appendix  B).  It  is  noted  that  for  all 
diluent  concentrations  Tav<  T^,  and  that  as  !Fay  approaches  the 
stability  limits  shift  toward  stoichiometric. 

In  Figure  20  similar  temperature  data  are  presented  for  CH^- 
dilute  oxygen  and  CH^-air.  The  T  given  for  CH^-air  is  the  temper¬ 
ature  at  stoichiometric.  For  CH^-(0.402+0.6N2)  and  a  given  chamber 
length  the  measured  combustion  temperatures  at  the  stability  limits  are 
approximately  the  same.  For  L  ■  J6  cm  Tfty  ■  2050i35°K>  and  for  L  ■ 

53  cm  Tay  -  2060t20°K.  For  CH^-(0.U02+0.6N2>  Tay  <*ad*  For  CHjj-air 
the  measured  combustion  temperature  at  stoichiometric  is  somewhat  less 
than  the  adiabatic  combustion  temperature,  but  approximately  the  same 
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as  Tav  for  ^-(O.liOg+O.dl^). 

To  facilitate  the  comparison  of  the  experimental  and  theoretical 
results  a  brief  summary  of  some  important  conclusions  obtained  from  the 
Sirignano-Crocco  analysis  is  mad*.  For  a  given  propellant  combination  the 
longitudinal  stability  limits  are  defined  by  f  -  Tcrlt,  where  T  ^  ie  the 
solution  of  equation  (3-4) .  Whenever  Tcrlt.<  Tm,v  two  regions  of  unstable 
combustion  will  exist,  one  on  either  3lde  of  stoichiometric.  Each  of 
these  regions  is  bounded  by  a  stability  limit  and  a  combustion  limit. 
Whenever  Tcrlt  fc  T[My  unstable  combustion  occurs  over  the  entire  range  of 
combustible  equivalence  ratios.  Thus  the  theory  provides  a  criterion  for 
the  number  and  location  of  unstable  regions. 

The  nearly  constant  value  for  T  (for  a  given  propellant 

av 

combination)  in  Figures  10  and  20  strongly  supports  the  concept  of  a 

critical  combustion  temperature  at  the  stability  limits.  The  plight 

dependence  of  T  on  L  is  discussed  in  Section  H.  The  theoretical 
av 

criterion  for  the  number  and  location  of  unstable  regions  is  also 
supported  by  experiment.  From  Figures  18  and  20  it  is  noted  that  when¬ 
ever  Tav  <  Tft(j  two  unstable  regions  exist,  one  on  either  side  of 
stoichiometric-  viz.  Hg-air,  Hg-dilute  air  and  CH^-Co.UOg-KJ^Ng) .  For 
CH^-air  the  measured  combustion  temperature  at  stoichiometric  is 

approximately  the  same  as  T  for  CH,  -(O.UO  +0.6N  ),  and  one  unstable 

av  4  2  2 
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reglon  (located  around  stoichiometric)  la  observed. 

The  upper  stability  Halts  of  the  oxygen-rich  unstable  regions 
for  H^-alr  and  H^-dllute  air  are  not  predicted  by  the  Slrlgnano-Crocco 
analysis.  The  probable  explanation  of  these  "additional"  stability 
limits  is  given  in  Section  D. 


Combining  the  stability  limit  data  of  Figures  17  and  19  vlth 

the  y  *8  for  the  tvo  propellant  combinations  it  is  possible  to  calculate 

(E/RT)llmit  from  equation  (3-4).  Using  these  values  for  (E/RT)linit  and 

the  values  for  Tay  from  Figures  18  and  20  it  is  possible  to  evaluate  Eiinif 

These  results  together  with  overall  activation  energies  for  the  tvo 

propellant  combinations  taken  from  the  chemical  kinetics  literature  are 

tabulated  in  Figure  21.  The  results  are  based  on  a  chamber  length  of 

53  cm.  It  is  noted  that  the  values  for  for  H2"air  and  air 

are  approximately  the  same.  This  result  supports  the  earlier  contention 

that  nitrogen  has  little  effect  on  the  activation  energy.  The  values  for 

E_.  obtained  from  (3-4)  are  approximately  half  the  corresponding 
limit 

activation  energies  obtained  from  the  literature. 


E 

E 


act  ^2^ 
limit  <CH4) 


Eact  Gty 


=  0.38 


*0.52 


However,  the  ratio  of  the  Ellalt  tor  the  two  reaction*  compares  favorably 
vith  the  ratio  of  the  activation  energies. 


E  limit 

E  limit  <»¥ 


*037 


E.ct  <h2> 
Eoet  <CH,» 


*  0.5t 


If  the  E  were  used  with  the  tabulated  values  of  (e/RT),,  .. 

act  limit 

to  calculate  f  ..  it  is  found  that  T  ..  >  T  for  both  propellant 
crit  crit  max 

combinations.  In  this  calculation  it  is  necessary  to  correct  the  tabulated 

(e/rT)..  for  the  temperature  dependence  of  Y  .  From  the  theory  when- 
limit 

ever  Tcr^  >  T,^  unstable  combustion  will  exist  over  the  entire  range 
of  combustible  mixture  ratios. 


There  are  several  possible  explanations  for  the  discrepancy 
between  the  magnitudes  of  Ei^m<f.  and  First  it  is  noted  that  the 

presence  of  stabilizing  effects  in  the  gas  rocket  (heat  transfer, 
friction,  finite  combustion  zone)  are  not  included  in  the  Sirignano- 
Crocco  analysis.  When  stabilizing  effects  are  present  the  perturbation 
in  the  energy  release  rate  required  to  support  an  oscillation  of  given 
strength  will  be  greater  than  when  no  stabilizing  effects  are  present. 
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Expressing  the  perturbation  in  the  energy  release  rate  in  terms  of  the 
pressure  perturbation  (assuming  an  lsentropic  wave) 

_l!  «  _e_  /id.\/PL\  (>*-2) 

f  rT  \  r  A  P I 

From  equation  (4-2)  It  Is  seen  that  for  a  given  pressure  perturbation  to 
achieve  an  increase  in  the  rate  perturbation  (for  a  given  propellant  com¬ 
bination)  a  decrease  in  7  is  required.  Thus  the  net  effect  of  stabilizing 
influences  is  a  shift  in  the  stability  limits  toward  lower  combustion 
temperatures.  In  later  sections  it  will  be  seen  that  stabilizing  effects 
are  not  sufficient  to  account  for  the  discrepancy  between  E.  and  E  . 


An  examination  of  the  combustion  zone  analysis  employed  in  the 

Sirignano-Crocco  theory  provides  a  possible  explanation  for  the  discrepancy 

between  E  and  E  .  In  the  Sirignano-Crocco  analysis  it  is  assumed 
limit  act 

that  the  combustion  process  in  the  gas  rocket  is  effectively  a  premixed 

turbulent  flame  in  which  the  local  energy  release  rate  is  given  by  equation 

(3-1).  The  total  energy  release  rate  (r  )  is  obtained  by  integrating 

total 

the  local  energy  release  rate  over  the  volume  of  the  combustion  zone.  For 
one-dimensional  flow 

s 

'total  *  /  rlx)Adx 

0 


where  r(x)  *  local  energy  release  rate,  8  *  combustion  zone  thickness  and 
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A  -  cross-cectional  area  of  the  combustion  chamber.  In  the  Sirignano-Crocco 
theory  it  is  assumed  that  8  is  independent  of  time  and  that  r(x)  can  be 
replaced  by  a  space-mean  energy  release  rate,  r  ,  which  is  a  function  of 
time.  Hence 

rl0t0l  ~  fAV  * 

Expressing  the  space-mean  energy  release  rate  in  terms  of  equation  (3-1), 


In  reality  the  combustion  zone  thickness  is  a  function  of  time  so  that 


r  ~  8(t)  exp 
total 


(4-4) 


For  a  premixed  turbulent  flame  the  combustion  zone  thickness  may  be  written 

(*5) 

8  « 

where  p  =  density  of  the  unburned  propellant,  Q  ■  eddy  diffusivity 

and  r  ■  space-mean  turbulent  reaction  rate.  The  eddy  diffusivity  may  be 
t 
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vritten 


$  ~uD 


where  u  ■  gae  velocity  tnd  D  ■  diaaeter  of  the  combustion  chamber.  The 
density  may  be  written 


P  ~ 


vhere  m  is  the  mass  injection  rate  (a  constant). 


Hence 


If  it  Is  assumed  that  the  turbulent  reaction  rate  can  be  expressed  as 


then  equation  (4-4)  may  be  written 

rtotol  ~  exp  I"  2  RTav}  • 

Following  the  earlier  approximation  of  Sirignano  and  Crocco  n  is  set  equal 


rtotol 


(4-5) 


to  zero 


expfik:-. 

AV  J 


For  the  energy  release  lav  M--5)  the  defining  equation  for  the  stability 
Halts  is 

(4-6) 

An  alternate  model  for  the  gas  rocket  combustion  process  has  been 
suggested  by  Slrlgnano.  In  this  model  it  Is  assumed  that  each  injection 
orifice  produces  a  conical-type  premixed  turbulent  flame  and  that  the 
combustion  zone  consists  of  a  ''sheet"  of  these  individual  flames.  For  this 
model  the  total  energy  release  rate  Is  given  by 


Sotal 


vhere  p  »  density  of  the  unburned  propellant,  q  ■  heat  of  reactlon/unlt 
mass,  =  turbulent  flame  speed  and  ■  area  of  the  flame  front.  The 
turbulent  flame  speed  is  a  function  of  time.  I*.  follows  that  the  flame  area 
is  a  function  of  time.  If,  as  a  first  approximation,  it  is  assumed  that  the 
density  and  velocity  of  the  approach  flow  are  constant  and  that  variations 
in  flame  area  (due  to  variations  in  flame  speed)  are  'negligible  then 


T 


r  ~  s 
total 


For  a  premixed  turbulent  flame  the  flame  speed  may  be  vrltten  (45) 


*T  ~  \/^T  • 

Proceeding  as  on  page  51  it  is  found  that 


r 


total 


(4-5) 


and  that 


(4-6) 


Thus  both  combustion  models  lead  to  similar  results.  Based  on  present 
knowledge  it  is  not  possible  to  determine  the  appropriate  model  for  the  gas 
rocket  combustion  process.  A  logical  extension  of  the  present  study  would 
include  experiments  designed  to  define  more  precisely  the  nature  of  the 
combustion  process  in  the  gas  rocket. 


The  E-j.imit  °t)taine<*  ^rom  equation  (4-6)  are  tabulated  in  Figure  21 
for  the  various  propellant  combinations.  It  is  seen  that  equation  (4-6) 
provides  E  ^  which  are  in  relatively  good  agreement  with  E  . 

2.  Dependence  of  Stability  Limits  on  £  (T  =  constant) 

Khskan  and  Browne  (37),  Friedman  and  Nugent  (38)  and  Lewis  and 
von  Elbe  (39)  have  reported  a  substantial  increase  in  the  rate  of  the  CO-O^ 
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reaction  by  the  addition  of  trace  amounts  of  H^.  This  reaction  was  employed 
in  an  experiment  designed  to  provide  information  on  the  dependence  of  the 
stability  limits  on  the  reaction  kinetics  (in  particular,  on  the  overall 
activation  energy).  In  the  experiment  the  regions  of  unstable  combustion 
for  CO-air  were  located  in  the  $-L  plane  with  and  without  the  presence  of 
the  catalyst.  The  shift  in  the  stability  limits  produced  by  the  addition 
of  the  Hg  was  correlated  with  the  Sirignano-Crocco  theory  and  available  rate 
data.  Air  was  used  as  the  oxidizer  primarily  because  it  was  found  that  for 
CO-Og  stable  combustion  existed  over  the  entire  range  of  combustible  mixture 
ratios-  ie.  the  stability  limits  fell  outside  the  combustion  limits.  Thus  it 
would  not  have  been  possible  to  obtain  the  required  data  using  CO-O^.  For 
CO-air  (with  correspondingly  lower  combustion  temperatures)  the  stability 
limits  shifted  toward  stoichiometric  and  fell  inside  the  combustion  limits. 

An  additional  reason  for  using  CO-air  was  the  severe  flashback  problem 
encountered  with  CO-Og. 

It  was  found  that  it  was  difficult  to  contain  the  reaction  inside 
the  combustion  chamber  when  using  dry  CO-air  (due  to  the  relatively  low 
flame  speed) .  To  decrease  the  probability  of  blow-off  the  flame  speed  was 
increased  by  the  addition  of  trace  amounts  of  water  vapor  to  the  CO  (4l). 

It  was  also  found  that  it  was  not  possible  'to  ignite  CO(HgO)-air 
mixtures  using  the  spark  ignition  system.  To  circumvent  this  problem  a  run 

procedure  was  employed  in  which  combustion  was  initiated  using  H  and  air. 

2 

Once  combustion  had  been  established  the  CO  and  air  flow  rates  were  set  to 
their  steady-state  values  (determined  by  mixture  ratio  and  combustion  pressure 


requirements).  The  It,  flow  used  to  initiate  combustion  was  set  so  that  when 
steady-state  operation  vas  obtained  the  flow  rate  was  that  required  to 
accelerate  the  reaction  (~  2.5#  molar  of  the  total  flow  rate).  Thus  the 
stability  of  the  combustion  with  and  without  the  was  determined  simply 
by  turning  the  flow  on  and  off.  By  repeating  this  procedure  for  various 
mixture  ratios  and  a  fixed  combustion  chamber  length  it  was  possible  to 
determine  the  shift  in  stability  limits  produced  by  the  addition  of  Hg. 

The  results  of  this  experiment  for  a  chamber  length  of  104  cm  and 

a  mean  combustion  pressure  of  J.8  atm  are  shown  in  Figure  22a.  The  stability 

limit  for  CO^Oj-air  is  <J>  =  2. 35*0.05*  This  corresponds  to  an  adia- 

d  limit 

batic  combustion  temperature  at  the  stability  limit  of  Tcrit  =  1173*20°K. 

The  stability  limit  for  CO^Oj-air-H^  is  $  ^  =  2-9010.05.  This 

—  o 

corresponds  to  a  T  .  =  1012113  K. 

crit 

As  a  first  step  in  the  interpretation  of  the  experimental  results 
the  effect  of  H^-addition  on  the  combustion  temperature  of  the  C0(H^0)-air 
reaction  was  determined.  As  noted  earlier  the  flow  rate  was  2. 5*1.0# 

(molar)  of  the  total  flow  rate.  Theoretical  calculations  indicate  that  this 
amount  of  will  increase  the  adiabatic  combustion  temperature  by  less  than 
12°K  over  the  range  of  equivalence  ratios  0.5  to  h.O.  This  conclusion  is 
supported  by  combustion  pressure  records  during  stable  combustion-  no  change 
in  P  was  observed  with  the  addition  of  the  Hg.  Both  the  increased  mass  flow 
and  the  increased  combustion  temperature  which  occur  with  the  addition  of 
serve  to  increase  the  combustion  pressure.  If  either  or  both  the  mass  addition 
or  temperature  effects  were  significant  a  change  in  P  would  be  observed. 
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Combining  the  stability  limit  data  of  Figure  22a  with  the  y  's 

of  the  combustion  products  it  is  possible  to  evaluate  E  for  CO(HoO)- 

limit 

air  and  CC^HgO^air-Hg.  The  obtained  from  combustion  laws  (3-l) 

and  (4-5)  and  the  overall  activation  energies  for  the  CC^HgOj-C^  reaction 
taken  from  the  literature  are  tabulated  in  Figure  22a.  As  was  true  for 
both  H2«air  and  CH^-Og  the  value  for  Ej.init  for  CO(HgO)-air  uslnfi 
combustion  law  ( 3—1)  is  approximately  half  the  activation  energy  obtained 
from  the  literature. 


E  limit  (CO) 
Eoct  (CO) 


0.56  . 


The  E  obtained  using  combustion  law  (4-5)  is  in  good  agreement  with 

limit 

the  literature  values  for  E  . .  From  Figure  22a  it  is  seen  that  E. . 

act  0  limit 

for  the  CC^HgOj-air  reaction  is  substantially  reduced  by  the  addition  of 
trace  amounts  of  H^. 


E  limit  (withH2> 
Elimit (wi,hout  H2> 


■  0-85  ±  0  02 


The  results  of  the  Sirignano-Crocco  analysis  will  now  be 
expressed  in  a  form  which  permits  comparison  with  kinetics  data.  The  mean 
energy  release  rate  (r)  can  be  expressed  in  two  ways  depending  on  the 
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combustion  lav  used  to  describe  the  gas  rocket  combustion  process- 
equation  (3“l)  or  (4-5).  In  terms  of  equation  (3“l) 


?  *  constant  exp 


and  in  terms  of  equation  (4-5) 


r  =  constant 


The  mean  energy  release  rate  is  related  to  the  mean  reaction  rate  (?)  by 


(4-7) 


Equation  (4-7)  can  be  used  to  compare  the  rates  of  the  C0(H2O)-air 
reaction  (subscript  l)  and  the  CC^I^Oj-air-Hg  reaction  (subscript  2), 


i  5  i  .  _9 2 

h  h  <»i 

Substituting  for  r  from  equation  ( 3— l) 


(4-8a) 
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and  substituting  for  r  from  equation  (4-5) 


(4-8b) 


To  determine  the  effect  of  addition  on  the  rate  of  the  CO(HgO)-air 

reaction  equations  (4-8)  aie  specialized  to  the  case  of  <X>  =  constant. 

In  view  of  the  earlier  discussion  of  the  effect  of  II  addition  on  the 

2 

combustion  temperature  of  the  C0(Hg0)-air  reaction  it  is  reasonable  to 
assume  (for  constant  )  that  q1  qg  and  Tg.  Hence  equations 

(4-8)  reduce  to 


(4-9a) 


(4-9b) 


In  equations  (4-9)  the  values  for  E,  and  E0  are  taken  to  be  the  E. . 

i  c  limit 

for  C0(Hg0)-air  and  C0(Hg0)-air-Hg  respectively.  Since  the  values  for 
E  in  equation  (4-9a)  are  half  those  for  E  in  equation  (4-9b)  it  is  seen 
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that  equations  (4— 9a)  and  (4-9b)  are  the  same.  Substituting  for  E^  E 
and  R  in  equations  (4-9) 

^p^sTwieoj  (4,0) 

By  inserting  a  value  for  T  in  equation  (4-10)  ?  /5^  is  determined.  The 
from  e<Juation  (4-10)  for  7  =  1732®K  is  tabulated  in  Figure  22b. 

It  is  emphasized  that  the  ?^/?2  obtained  from  (4-10)  is  the  ratio  of  the 
mean  reaction  rates  derived  from  the  Sirignano-Crocco  analysis  and  the 
experimentally  observed  shift  in  the  stability  limits. 

It  is  of  interest  to  compare  the  ?  /?  obtained  from  equation 
(4-10)  with  available  kinetics  data  for  the  CO-air  reaction.  Two  sources 
of  kinetics  data  are  the  laminar  flame  speed  data  of  Jahn  (39)  and  the 
laminar  flame  structure  studies  of  Friedman  and  Nugent  (38).  The  flow 
in  the  gas  rocket  is  highly  turbulent  hence  the  combustion  process  in  the 
gas  rocket  is  essentially  a  premixed  turbulent  flame.  To  compare  the 
results  of  equation  (4-10)  with  the  kinetics  data  (38,  39)  it  is  necessary 
to  infer  turbulent  reaction  rates  from  laminar  flame  data. 


\ 


There  are  two  basically  different  models  for  premixed  turbulent 


flames.  The  first  of  these,  termed  the  wrinkled  flame  model,  assumes 
that  a  turbulent  flame  is  simply  a  laminar  flame  front  which  has  been 

"wrinkled"  due  to  local  variations  in  the  flow  velocity  (42,  43) .  The 
wrinkled  flame  model  is  generally  applied  to  flows  with  low  intensity, 
large  scale  turbulence.  The  second  model,  termed  the  distributed 
reaction  zone  model,  describes  a  turbulent  flame  as  a  "zone  of  reaction 
distributed  in  depth,"  making  no  direct  reference  to  laminar  flames  (44). 
The  distributed  reaction  zone  theory  is  generally  applied  to  flows  with 
high  intensity,  small  scale  turbulence. 

Considering  the  intensity  and  scale  of  the  turbulence  in  the 
gas  rocket  it  is  felt  that  the  distributed  reaction  zone  model  best 
describes  the  gas  rocket  combustion  process.  The  distributed  reaction 
zone  model  assumes  that  the  propagation  of  a  premixed  turbulent  flame 
is  analogous  to  the  propagation  of  a  premixed  laminar  flame  (44,  45). 
Hence 


(4-11) 


The  differential  equation  describing  the  turbulent  flame  is  similar  to 
the  equation  describing  a  laminar  flame,  with  the  laminar  transport 
properties  being  replaced  by  turbulent  transport  properties.  The  solution 
of  this  differential  equation  provides  a  relationship  between  the  turbulent 
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flame  speed  (S  )  and  the  laminar  flame  speed  (S  ), 
t  1 

Vt  5l*l 

$  *  a 


(4-12) 


where  a  =  thermal  diffusivity,  3  -  thickness  of  the  turbulent  flame, 

t 

and  3  *  thickness  of  the  laminar  flame.  Combining  equations  (4-ll) 

and  (4-12)  an  expression  for  the  turbulent  reaction  rate  is  obtained. 


(4-13) 


From  the  theory  of  the  propagation  of  premixed  laminar  flames  (45) 


(4-14) 


where  =  mean  laminar  reaction  rate.  Combining  equations  (4-13)  and 
(4-14), 


(^-15) 


! 


Equation  (4-13)  or  (4-15)  can  be  used  with  the  kinetics  data 
(38,  39)  to  obtain  a  ratio  of  reaction  rates  for  the  CO-air  (subscript  l) 
and  CO-air-H^  (subscript  2)  reactions. 


fikif 

^  \§L2  / 

2 

(£l  (1l\ 

\W\ai)  v?L2  I 


(4-16) 


The  addition  of  trace  amounts  of  to  the  CO-air  reaction  produces  only- 
slight  variations  in  the  combustion  temperature  and  in  the  density  and 
composition  of  the  reactants  end  products.  Thus  the  addition  should 
have  only  a  slight  effect  on  the  eddy  and  thermal  diffusivities-  ie. 


t 

r 


S|  »  S; 


fl|  «  a2 
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Slnce  the  addition  of  trace  amounts  of  H2  accelerates  the  CO-air  reaction 
it  Is  to  be  expected  that  the  thickness  of  both  the  laminar  flame  and  the 
turbulent  flame  will  decrease.  It  is  felt,  however,  that  the  ratio 

8^/  8  will  be  only  slightly  affected  by  the  addition  of  Hg.  Hence, 
as  a  first  approximation  equation  (4-l6)  can  be  written 


Jahn  (39)  has  measured  the  effect  of  addition  on  the  laminar 
flame  speed  in  premixed  CO(HgO)-air  at  atmospheric  pressure.  Jahn's  data 
are  used  with  equation  (4-17a)  to  estimate  r^/?^  ?or  the  CO(B^O)-air 
reaction  at  a  temperature  of  1732°K  (  $  ■  1.49).  The  result  is  tabulated 

in  Figure  22b. 

Friedman  and  Nugent  (38)  have  studied  the  effect  of  the  addition 
of  trace  amounts  of  Hg  and  water  vapor  on  the  CO-Og  reaction.  Overall 
specific  rate  constants  were  inferred  from  temperature  and  concentration 
traverses  through  premixed  laminar  flames.  They  also  measured  the  effect 
of  Hg  and  HgO  addition  on  the  laminar  flame  speed.  The  studies  were 
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performed  for  a  pressure  of  30  mm  Hg.  Friedman  and  Bugent  found  that 
the  C0-02  reaction  was  first  order  with  respect  to  the  concentration 
of  CO. 

£fei«-K[CCj  (4-lfl) 

The  factor  k  in  equation  (4-l8)  is  the  overall  specific  rate  constant 
measured  by  Friedman  and  BUgent.  Using  equation  (4-18)  and  the  rate 
constant  data  (referred  to  a  temperature  of  1732*K)  is  calculated. 

This  value  of  is  substituted  into  equation  (4 -17b)  to  obtain 

?ti/?t2«  Since  the  value  of  the  specific  rate  constant  was  different  as 
determined  from  temperature  traverse  data  and  composition  traverse  data 
two  values  of  are  tabulated  in  Figure  22b.  Using  equation  (4-17a) 

and  the  flame  speed  data  (referred  to  a  temperature  of  1732*K)  7^/7^ 
was  calculated.  This  7^/7^  1b  tabulated  in  Figure  22b. 

The  data  of  Figure  22b  allow  a  comparison  to  be  made  between  the 
rate  ratio  obtained  from  the  Sirlgnano-Crocco  theory  and  instability  data 
and  the  ratio  obtained  from  kinetics  data  and  a  model  for  the  gas  rocket 
combustion  process.  In  making  this  comparison  it  is  important  to  under¬ 
stand  some  of  the  limitations  of  the  results  of  Figure  22b.  The  equations 
used  to  relate  the  turbulent  reaction  rate  to  the  laminar  reaction  rate 
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and  the  laminar  flame  speed  are  order  of  magnitude  approximations.  In 
addition  the  studies  of  Friedman  and  Nugent  (3d)  have  shown  that  the 
laminar  flame  speed  and  rate  constant  depend  on  the  concentration  of 
water  vapor.  Since  the  water  vapor  content  In  the  Instability  experi¬ 
ments  and  In  the  various  kinetics  experiments  (38,  39)  vere  not  the 
same  differences  In  the  reaction  rates  are  to  be  expected.  Furthermore, 
the  data  of  References  38  and  39  were  obtained  at  pressures  which  were 
significantly  smaller  than  the  combustion  pressure  In  the  gas  rocket, 
ftiere  are  no  data  available  on  the  effect  of  pressure  on  the  relatively 
complex  CO-O^  reaction.  Water  vapor  and  pressure  effects  should  not 
have  a  strong  effect  on  the  results  of  Figure  22b  due  to  the  fact  that 
ratios  of  reaction  rates  sure  used  (and  such  effects  would  tend  to  cancel). 
Considering  these  and  other  limitations  it  la  possible  that  the  relatively 
good  agreement  between  the  results  of  the  Sirignano-Crocco  theory  and 
kinetics  data  is  fortuitous.  However,  several  Import  sunt  conclusions  can 
be  drawn  from  the  data  of  Figure  22b.  The  Sirignano-Crocco  theory  (using 
instability  data)  and  chemlcsd  kinetics  studies  both  show  an  increase  in 
the  rate  of  the  C0-02  reaction  with  the  addition  of  trace  amounts  of  H^. 
The  magnitude  of  the  rate  increase  obtained  from  the  kinetics  studies 
is  of  the  same  order  as  that  given  by  the  Sirignano-Crocco  theory  using 
instability  data. 

In  summsu’y  it  is  emphasized  that  the  experiments  described  in 
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thls  section  afford  a  critical  test  of  a  chemical  kinetic  Instability 
driving  mechanism.  The  two  important  kinetic  parameters,  temperature 
and  activation  energy,  have  been  varied  independently  :Ln  instability 
experiments,  and  the  observed  changes  in  the  stability  limits  have  been 
shown  to  correlate  well  with  am  instability  theory  based  on  a  chemical 
kinetics  driving  mechanism. 

D.  Dependence  of  Instability  Strength  on  T 

For  a  given  M  the  Sirignano-Crocco  analysis  gives  the  strength 
of  the  instability  (  &p/p)  as  a  function  of  E/RT  and  Y .  For  a  given 

propellant  combination  E  is  approximately  constant  and  Y  is  a  function 

* 

only  of  T.  Hence  for  given  propellants  A  P/P  is  a  function  only  of  T. 

To  determine  the  temperature  dependence  of  A  ?/V  an  experiment 
was  carried  out  in  which  AP/F  was  measured  as  a  function  of  $  for  a 
given  propellant  combination  (holding  L,  P  and  R  constant).  For  a  given 
combustion  pressure  T  is  a  function  of  $  only,  hence  the  experiment 
provides  A  P/F  as  a  function  of  T, 

The  results  of  a  typical  experiment  using  H^-air  ere  presented 
in  Figure  23  (in  which  Ap/F  is  plotted  as  a  function  of  the  mean 
adiabatic  combustion  temperature).  The  data  were  obtained  for  L  -  38  cm, 
P  ■  7.0±O:3  atm,  R  -  0.015±0.002,  the  31_bole  showerhead  injector  and 


*  Over  the  range  of  pressures  of  interest  Y  is  independent  of  pressure. 
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the  oxidizer-rich  unstable  region.  The  fuel -rich  unstable  region  vas  too 
narrow  to  permit  satisfactory  measurement  of  shock  strength.  The  mean 
Mach  number  of  the  combustion  gases  vas  calculated,  using  known  mass  flow 
rates  and  average  values  for  the  combustion  temperatures.  The  heavy 
vertical  bar  gives  the  variation  in  shock  strength  during  a  given  run 
(due  to  non-periodic  fluctuations). 

From  Figure  23  it  is  seen  that  as  T  decreases  &P/P  Increases 
to  a  maximum  and  then  decreases.  For  sufficiently  low  T  the  shock 
strength  decreases  to  zero  producing  the  upper  stability  limits  observed 
for  the  oxidizer-rich  unstable  regions  for  H^-air  and  H^-dilute  air 
mixtures  (Figure  IT).  The  fact  that  the  upper  stability  limit  is  not 
clearly  defined  is  due  to  the  relatively  slow  decrease  in  shock  strength 
with  T.  For  relatively  weak  shock  waves  it  is  difficult  to  distinguish 
the  combustion  oscillations  from  the  random  fluctuations  of  turbulence. 
Hence  when  the  shock  strength  decreases  slowly  precise  location  of  the 
stability  limits  is  not  possible. 

The  results  of  an  experiment  performed  by  Pelmas  (27)  using 
CHj,-(0.4c>2+0.6N2)  are  presented  in  Figure  24.  These  data  were  obtained 
for  L  ■  79  cm,  V  ■  7.8  atm,  a  porous  plug  injector  and  the  fuel-rich 
unstable  region.  Since  H  was  not  specified  in  Reference  27  its  value 
was  calculated  using  estimated  flow  rates,  the  known  P  and  the  adiabatic 
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combustion  temperature.  From  this  calculation  M  was  estimated  to  be 
0.015.  The  results  of  Figure  2k  agree  qualitatively  with  the  results 
of  Figure  23  in  that  as  T  decreases  A  ?/$  increases  to  a  maximum  and 
then  decreases. 

In  contrast  to  the  experimental  results  the  Sirignano-Crocco 
analysis  predicts  that  Ap/p  increases  monotonically  with  decreasing 
T  (See  Figure  14).  Hie  difference  between  experimental  and  analytical 
results  can  be  explained  as  follows.  As  the  strength  of  the  instability 
Increases  the  losses  in  the  wave  increase,  and  hence  the  energy  supplied 
to  the  oscillation  by  the  combustion  process  must  increase.  At  lower 
combustion  temperatures  the  mass  fraction  of  combustibles  in  the  react¬ 
ant  mixture  decreases,  and  hence  the  energy  which  can  be  supplied  by 
the  combustion  process  (for  constant  mass  injection  rate)  decreases. 

It  is  argued  that  below  some  T  the  energy  supplied  by  the  combustion 
process  is  not  sufficient  to  support  an  oscillation  of  the  strength 
required  by  the  theory.  Further  since  the  energy  supplied  by  combus¬ 
tion  decreases  with  the  combustion  temperature  it  follows  that  the 
strength  of  the  instability  must  decrease  as  T  decreases  below  the  above 
value. 


The  theoretical  variation  of  CkT/V  with  T  (based  on  combus¬ 
tion  law  3-l)  Is  plotted  in  Figure  23  for  H^-air  with  M  =  0.015  and 


E  «*  15*1  kcal/mole  (the  value  for  E  obtained  from  experiment)  and 

limit 

37.8  kcal/mole  (the  activation  energy).  A  similar  curve  is  plotted  in 
Figure  24  for  CH^-(0.402+0.6N2)  with  M  =  0.015  and  E  =  E^  =  23.6 
kcal/mole.  The  theoretical  curves  based  on  combustion  law  (4-5)  corre¬ 
spond  to  the  E  =  15.1  kcal/mole  curve  (Figure  23)  and  the  E  =  23.6 
kcal/mole  curve  (Figure  24).  The  slopes  of  the  theoretical  curves  are 
significantly  smaller  than  the  initial  slopes  of  the  experimental  curves. 
From  Figure  23  it  is  seen  that  this  result  cannot  be  attributed  to  un¬ 
certainties  in  the  value  for  E.  Furthermore,  the  slight  uncertainties 
in  M  are  not  sufficient  to  account  for  the  discrepancies  between  exper¬ 
imental  and  theoretical  results.  It  is  probable  that  the  discrepancies 
are  the  result  of  the  simplifications  inherent  in  the  combustion  zone 
analysis  and  energy  release  law  (See  Chapter  3). 

E.  Dependence  of  Instability  Strength  on  R 

For  a  given  propellant  combination  and  combustion  temperature 
the  Sirignano-Crocco  analysis  gives  the  instability  strength  as  a 
function  of  the  mean  Mach  number  of  the  combustion  gases.  To  determine 
the  dependence  of  the  instability  strength  on  M  an  experiment  was 
carried  out  in  which  A  P/P  was  measured  as  a  function  of  <D  (ie.  T) 
for  two  values  of  R  (for  given  L,  P  and  propellant  combination).  The 
variation  in  R  was  obtained  by  using  exhaust  nozzles  with  different 
flow  areas .  The  dependence  of  &  P/P  on  M  was  obtained  by  comparing 


the  instability  strengths  for  the  two  R  at  a  given  combustion  temperature. 

The  results  of  a  typical  experiment  using  Hg-air  are  presented 
in  Figure  25.  In  this  figure  &P/P  is  plotted  against  the  adiabatic 
combustion  temperature  for  M  =  0.015  and  0.030  (i0.002).  The  data  were 
obtained  for  L  s  38  cm,  7  ■  4.4±0. 3  atm,  the  31-hole  6howerhead  in¬ 
jector  and  the  oxidizer-rich  unstable  region.  The  curve  for  M  *  0.015 
was  obtained  using  the  converging-diverging  nozzle  with  ;  =  k6.  The 
curve  for  M  =  0.030  was  obtained  using  the  converging-diverging  nozzle 
with  CR  =  23.  The  M  were  calculated  using  the  known  mass  flow  rates  and 
average  values  for  the  combustion  temperatures.  The  heavy  vertical  bars 
indicate  a  non-periodic  variation  in  the  shock  strength  during  a  given 
run.  The  variation  in  shock  strength  is  larger  for  the  larger  value  of 
R.  This  result  is  to  be  expected  since  the  intensity  of  the  turbulent 
fluctuations  increases  as  the  mean  flow  velocity  (ie.  R)  increases. 
Qualitatively  the  Ap/F  curves  are  the  same  as  in  Figures  23  and  2h- 
ie.  with  decreasing  T  the  instability  strength  increases  to  a  maximum 
and  then  decreases.  The  shift  in  the  stability  limit  for  the  Cr  =  23 
nozzle  is  discussed  in  Section  J. 

In  Figure  2 6  the  experimental  and  theoretical  shock  strength 


ratio 
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is  plotted  as  a  function  of  the  adiabatic  combustion  temperature.  The 
experimental  shock  strength  ratio  is  obtained  from  the  data  of  Figure 
25.  The  theoretical  shock  strength  ratio  is  obtained  from  the 
Sirignano-Crocco  analysis  using  E  =  E^^  and  the  / 's  for  the  equilib¬ 
rium  composition  of  combustion  products.  Since  the  stability  limits  are 
different  for  the  tvo  M  the  are  different  for  the  two  M  : 

E.  .  4.  s  «  *  l1*-*2  kcal/mole  and  E,  .  JA|  -  .  =  13.6  kcal/mole. 

The  results  of  Figure  2 6  are  also  applicable  for  the  combustion  law  (U-5). 
Considering  the  uncertainties  in  the  instability  strength  it  is  seen 
that  there  is  good  agreement  between  theory  and  experiment  even  for  the 
range  of  T  where  the  shock  strength  decreases  with  the  combustion  tem¬ 
perature.  This  result  is  to  be  expected  since  M  enters  into  the  analysis 
independent  of  the  combustion  zone  analysis  or  the  form  of  the  energy 
release  law. 

F.  Harmonic  Mode  and  Transverse  Mode  Instabilities 

In  the  course  of  the  experimentation  it  was  found  that  longi¬ 
tudinal  combustion  oscillations  could  occur  in  a  form  other  than  a 


single  shock  wave  propagating  In  the  combustion  chamber.  Two  other 
forms  of  longitudinal  oscillations  have  been  observed-  one  in  which 
two  shock  waves  simultaneously  propagate  In  the  chamber  and  another  In 
which  three  shock  waves  simultaneously  propagate  In  the  chamber. 
Pressure-time  histories  and  x-t  diagrams  for  the  three  forms  of  longi¬ 
tudinal  oscillations  are  presented  In  Figures  27,  28  and  29.  The 
pressure-time  histories  were  obtained  at  a  station  approximately  16  cm 
downstream  from  the  injector  using  Hg-air  with  L  «  122  cm  and  P  *  7.8 
atm.  Figure  27  is  the  pressure-time  history  and  x-t  diagram  for  a 
single  shock  wave  oscillating  with  a  frequency  of  260  cps  (  $  =  3.25). 

Figure  28  is  the  pressure-time  history  and  x-t  diagram  for  two  shock 
waves  (  $  2.75).  Note  that  the  two  waves  travel  in  opposite 

directions,  one  wave  reflecting  from  the  nozzle,  the  other  from  the 
injector.  Hie  apparent  oscillation  frequency  is  520  cps.  Figure  29 
is  the  pressure-time  history  and  x-t  diagram  for  three  shock  waves 
(  <t>  ^  2.50).  Note  that  the  pressure-time  history  shows  a  combination 
of  the  two  shock  wave  and  three  shock  wave  instability  forms.  Hie 
apparent  oscillation  frequency  for  three  shock  waves  is  750  cps.  From 
Figures  27  through  29  it  is  seen  that  the  apparent  oscillation  frequen¬ 
cies  for  the  various  forms  of  longitudinal  instabilities  are  approxi¬ 
mately  integral  multiples  of  the  oscillation  frequency  of  the  single 


shock  wave. 
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where  fft  «  apparent  oec illation  frequency  of  the  n  mode,  f  = 
oscillation  frequency  of  a  single  shock  wave  and  n  =  an  Integer  =  1,  2 
or  3.  Bie  various  forms  of  longitudinal  oscillations  will  be  classified 
according  to  frequency  and  termed  the  fundamental  mode  (n  »  l),  second 
harmonic  (n  =  2)  and  third  harmonic  (n  ■  3)  by  analogy  with  acoustic 
oscillations.  Harmonic  mode  longitudinal  oscillations  of  the  type  shown 
in  Figures  27  through  29  have  been  observed  in  gas  rockets  by  several 
investigators.  Pelmas  (27)  reported  the  existence  of  harmonics  in  the 
vicinity  of  the  stability  limits  for  H^-air  and  CH^-Co.UC^+O.Sl^) . 

Tsuji  and  Takeno  (34)  reported  the  existence  of  harmonics  for  natural 
gas-  air. 


In  Figure  30  the  regions  of  the  various  harmonic  mode  insta¬ 
bilities  are  located  in  the  $  -L  plane  for  H^-air,  P  =  7*8  atm,  the 
31-hole  showerhead  injector  and  the  Cr  =  46  converging-diverging  nozzle. 
In  the  region  designated  "fundamental"  only  the  fundamental  mode  oscil¬ 
lation  (n  *  l)  occurs.  In  the  region  designated  "second  harmonic"  only 
the  second  harmonic  mode  (n  =  2)  oscillation  occurs.  There  is  some  over¬ 
lap  of  these  two  regions  (solid  cross-hatching),  and  in  this  overlap 
region  both  the  fundamental  and  second  harmonic  modes  are  observed. 
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For  a  given  equivalence  ratio  there  is  a  maximum  chamber  length  for  the 
appearance  of  the  fundamental  mode-  ie.  above  a  certain  chamber  length 
the  combustion  is  stable  with  respect  to  the  fundamental  mode.  The 
appearance  of  the  third  harmonic  mode  in  the  vicinity  of  the  stability 
limit  for  L  >  100  cm  indicates  that  there  is  probably  an  upper  length 
limit  for  the  appearance  of  the  second  harmonic  as  well.  Furthermore 
for  a  given  equivalence  ratio  the  order  of  the  harmonic  modes  (ie.  n) 
increases  as  the  chamber  length  increases-  ie.  longer  chambers  favor  the 
higher  harmonics.  For  a  given  chamber  length  the  order  of  the  harmonics 
increases  as  the  equivalence  ratio  approaches  stoichiometric-  ie.  the 
higher  harmorV.s  lie  closer  to  stoichiometric.  The  observations  concern¬ 
ing  the  ordering  of  the  harmonics  are  in  agreement  with  the  results  of 
Pelmas  (27)  and  Tsuji  and  Takeno  (3M» 

In  Figure  33  the  regions  of  the  various  harmonic  mode  instabil¬ 
ities  are  located  in  the  $  -P  plane  for  H^-air,  L  o  1 04  cm,  the  31-hole 
showerhead  injector  and  the  Cp  «  h6  converging-diverging  nozzle.  As  in 
Figure  3°  there  is  a  region  of  fundamental  mode  oscillations  and  a  region  of 
second  harmonic  mode  oscillations  (with  some  overlap  of  the  two  regions). 
For  a  given  equivalence  ratio  there  is  a  maximum  P  for  the  appearance  of 
the  fundamental  mode-  ie.  above  a  certain  mean  combustion  pressure  the 
combustion  is  stable  with  respect  to  the  fundamental  mode.  The  appearance 
of  the  third  harmonic  mode  in  the  vicinity  of  the  stability  limit  for 
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P  >  7«8  atm  indicates  that  there  is  probably  an  upper  pressure  limit  for 
the  appearance  of  the  second  harmonic  mode  as  well.  In  addition  for  a 
given  equivalence  ratio  the  order  of  the  harmonic  mode  increases  as  P 
increases-  ie.  higher  pressures  favor  higher  harmonics.  As  before  the 
order  of  the  harmonics  Increases  as  the  equivalence  ratio  approaches 
stoichiometric.  These  observations  are  in  agreement  with  the  results  of 
Tsuji  and  Takeno  (3*0  • 

The  Sirignano-Crocco  analysis,  by  assumption,  considers  only  the 
fundamental  mode  of  longitudinal  instability.  In  principle  their  analysis 
could  be  modified  to  include  harmonic  modes.  This  modification  would 
involve  altering  Assumption  2  (Chapter  3)  to  allow  for  more  than  one  shock 
wave,  and  obtaining  a  solution  subject  to  this  new  assumption.  Such  a 
modification,  if  possible,  would  provide  little  additional  information 
about  the  instability  phenomenon.  In  particular  since  the  mode  of  the 
oscillation  enters  the  analysis  by  assumption  it  would  not  be  possible 
to  obtain  a  criterion  for  the  occurence  of  the  harmonic  modes. 

A  qualitative  criterion  for  the  occurence  of  harmonic  mode 
instabilities,  based  on  the  fact  that  in  the  gas  rocket  the  combustion 
zone  has  a  finite  thickness,  has  been  devised.  For. a  combustion  zone 
of  finite  thickness  there  is  a  finite  time  between  the  injection  of  a 
gas  particle  and  its  final  conversion  to  burned  products.  If  T  is  the 
residence  time  of  a  gas  particle  in  the  combustion  zone  (overall  reaction 


time)  and  0  is  the  apparent  period  of  the  oscillation  then  the  criterion 
for  the  occurence  of  harmonic  mode  instabilities  can  be  expressed  as 


c  <  e/r  <  d 


(4-19) 


where  c  and  d  are  some  critical  parameters  (not  necessarily  constants). 
The  favored  harmonic  mode  is  the  one  for  which  0  /T  is  in  (c,d). 
Whenever  0  /X  exceeds  d  the  next  higher  harmonic  mode  (smaller  0  ) 
is  favored.  Whenever  0  /T  becomes  less  than  c  the  next  lower  harmonic 
mode  (larger  0  )  is  favored.  The  concept  of  a  critical  range  of  values 
for  0  /X  Js  supported  by  the  results  of  Tsuji  and  Takeno  (3*0-  for  a 
given  P  and  $  they  found  that  the  boundaries  between  regions  of  harmonic 
mode  oscillations  were  defined  by  0  «  constant.  The  criterion  (4-19) 
can  be  applied  to  explain  the  dependence  of  the  boundaries  separating 
the  various  harmonic  regions  on  L,  T  and  $  .  For  given  P  and  $  as  L 
is  increased  0  increases  and  T  remains  approximately  constant.  Hence 
as  L  increases  0  /T  increases.  If  the  increase  in  L  is  sufficient  to 
cause  0  / T  to  exceed  d  the  next  higher  harmonic  mode  is  fa\ored. 

Hence  (4-19)  predicts  that  the  higher  hirnnics  are  favored  at  longer 
chamber  lengths.  This  result  is  in  agreement  with  the  results  of  Figure 
3C.  For  given  L  and  $  as  P  is  increased  0  remains  approximately 
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constant  and  T  decreases  (reaction  rate  increases).  Hence  as  P  increases 
0  /T  increases.  If  the  increase  in  P  is  sufficient  to  cause  0  /  T 
to  exceed  d  the  next  higher  harmonic  mode  is  favored.  Hence  (4-19) 
predicts  that  the  higher  harmonics  are  favored  at  higher  combustion 
pressures.  This  result  ic  in  agreement  with  the  results  of  Figure  33* 

For  given  L  and  P  as  $  approaches  stoichiometric  (ie.  as  T  increases) 
both  and  T  decrease  (both  the  sound  speed  and  reaction  rate  increase 
with  increasing  T) .  Since  the  reaction  rate  is  a  stronger  function  of  T 
than  is  the  sound  speed  1*  decreases  faster  than  0  as  T  increases. 

Hence  as  $  approaches  stoichiometric  0  /T  increases.  If  0/1* 
exceeds  d  the  next  higher  harmonic  mode  is  favored.  Hence  (4-19)  predicts 
that  the  higher  harmonics  are  favored  closer  to  stoichiometric.  This 
result  is  in  agreement  with  the  data  of  Figures  30  and  33  and  with  the 
results  of  Pelmas  and  Tsuji  and  Takeno. 

The  harmonic  mode  criterion  given  by  equation  (4-19)  is  some¬ 
what  over-simplified.  However,  the  concept  of  a  critical  range  of  values 
for  the  parameter  0  /T  does  account  for  the  essential  features  of  the 
experimental  data,  and  indicates  that  to  account  for  harmonic  mode  in¬ 
stabilities  the  existence  of  a  finite  combustion  zone  mu6t  be  considered. 
Crocco  (33)  has  recently  advanced  a  theory  for  longitudinal  gas-phase 
instability,  assuming  a  combustion  zone  of  finite  thickness  and  employing 
a  rate  function  which  is  somewhat  more  realistic  than  ( 3- 1 ) • 


On  occasion  the  first  tangential  transverse  mode  would  occur 


simultaneously  with  the  longitudinal  mode  in  the  vicinity  cf  the  stabil¬ 
ity  limits.  A  typical  pressure-time  history  is  shown  in  Figure  31*  The 
frequency  of  the  longitudinal  mode  is  630  cps.  The  frequency  of  the 
first  tangential  mode  is  approximately  7500  cps. 

G.  Effect  of  Mean  Combustion  Pressure  on  the  Instability 

1.  Dependence  of  Stability  Limits  on  P 

To  determine  the  dependence  of  the  stability  limits  on  the  mean 
combustion  pressure  (P)  the  regions  of  unstable  combustion  for  H^-air  were 
located  in  the  $  -P  plane.  The  results  for  combustion  chamber  lengths  of 
33  cm  and  104  cm  are  shown  in  Figures  32  and  33  respectively.  These  results 
were  obtained  using  the  31-hole  showerhead  injector  and  =  46  converging- 
diverging  nozzle.  Only  the  oxidizer-rich  unstable  regions  are  shown  since 
variations  in  the  fuel-rich  stability  limits  over  the  range  of  P  inves¬ 
tigated  were  negligible.  For  the  range  of  $  investigated  (2.0  to  4.0) 
no  upper  stability  limit  was  observed  for  L  =  104  cm.  It  was  not  possible 
to  obtain  $>4.0  using  available  critical  flow  orifices.  Considering  the 
results  of  Figures  17  and  32  it  is  probable  that  the  upper  stability  limit 
exists  for  some  $  >4.0.  The  present  study  involves  only  the  limits  sepa¬ 
rating  regions  of  stable  combustion  from  regions  of  unstable  combustion. 

From  Figures  32  and  33  lb  is  seen  that  as  P  is  decreased  the  $  -range  for 
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unstable  combustion  decreases  and  the  values  for  $  ,  J  (at  the  lower 

limit 

stability  limit)  shift  away  from  stoichiometric.  Although  unstable 
combustion  was  encountered  over  the  entire  range  of  P  investigated,  the 
data  seem  to  indicate  that  there  is  some  minimum  value  for  P  below  which 
the  combustion  will  be  stable.  These  conclusions  are  in  qualitative  agree¬ 
ment  with  the  results  of  other  experimental  studies  of  combustion  insta¬ 
bility  in  gas  rockets  (31,  3*0*  In  these  studies  a  decrease  in  the  $  - 
extent  of  the  unstable  regions  with  decreasing  P  and  a  minimum  P  for 
unstable  combustion  were  observed.  Some  typical  results  from  Reference 
31  for  methane-air  are  shown  in  Figure  31** 

In  the  Sirignano-Crocco  analysis  the  stability  limits  are  in¬ 
dependent  of  P.  To  explain  the  differences  between  the  experimental 
observations  and  the  theoretical  results  the  effect  of  P  on  various 
combustion  parameters  will  be  examined,  and  the  dependence  of  the 
stability  limits  on  these  parameters  will  be  noted. 

Effect  of  P  on  the  Adiabatic  Combustion  Temperature 

As  P  decreases  the  overall  range  of  adiabatic  combustion 
temperatures  decreases.  This  effect  is  due  primarily  to  increased  dis¬ 
sociation  at  the  lower  pressures.  The  maximum  adiabatic  combustion 
temperatures  for  H^-air  for  various  pressures  are  shown  below. 
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As  noted  earlier,  a  decrease  in  the  overall  range  of  combustion  temper¬ 
atures  results  in  a  shift  in  the  stability  limits  toward  stoichiometric. 

Thus  a  decrease  in  P  should  produce  a  shift  in  $  ,  ,  ,  toward  stoi- 

]  imit 

chiometric.  This  effect  is  opposite  to  the  experimental  observations  of 
Figures  32  and  33.  It  is  concluded,  therefore,  that  the  dependence  of 
T  on  P  cannot  explain  the  observed  shift  in  stability  limits. 

Effect  of  P  on  y 

Equation  (3-4)  defines  the  stability  limits  in  terms  of  7  . 

If  7  is  a  function  of  pressure  then  (E/FT)  ,  and  hence  the  stability 

limit 

limits,  is  a  function  of  pressure.  7  has  been  computed  for  the  equilib¬ 
rium  products  of  combustion  of  H^-air  for  various  P  and  $  .  Over  the 

range  of  P  of  interest  and  for  $  in  the  vicinity  of  4  V 

1  imit 

assumes  a  constant  value  of  1.29«  It  is  concluded,  therefore,  that  the 
shift  in  stability  limits  with  P  cannot  be  explained  in  terms  of  the 


pressure  dependence  of  Y  . 


Effect  of  P  on  Heat  Transfer 

The  rate  of  heat  transfer/unit  area  to  the  walls  of  the  combus¬ 
tion  chamber  (Q)  may  be  approximated  by 

£  »hTg 

where  h  =  convective  heat  transfer  coefficient  (defined  by  equation  4-20) 
and  T  =  mean  gas  temperature.  The  heat  transfer  coefficient  can  be 
expressed  as  (46) 


h_D  s 

K 


constont 


/  .  \0-8 

te) 


'l±Coj 


(4-20) 


where  D  =  diameter  of  the  chamber,  K  -  thermal  conductivity  of  the 
combustion  gases,  m  =  mass  flow  rate,  jj,  =  viscosity  of  the  gases  and 
Cp  =  specific  heat  at  constant  pressure  of  the  gases.  Neglecting  the 
pressure  and  temperature  dependence  of 


H- 


,  cp  and  /C  , 


h  »  constant  m 


0.8 


Hence , 
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(£  »  constont  m  Tq. 


For  critical  flow  at  the  exhaust  nozzle 


rh  =  constont 


Hence , 


yTg  . 


(£  «  constont  P0^  Tg0  6  . 


(4-21) 


From  equation  (4-21)  It  is  seen  that  as  P  increases  Q  increases-  ie.  heat 
losses  increase.  From  equation  (4-3)  it  is  seen  that  as  heat  losses  in¬ 
crease  the  stability  limits  tend  to  shift  to  lower  combustion  temperatures- 
ie.  away  from  stoichiometric.  The  decrease  in  gas  temperature  due  to  heat 
transfer  is  given  by 


Tod  -  Tg  = 


m£p  . 


Neglecting  the  temperature  dependence  of  c^  this  equation  may  be  written 


Tod  -  Tg  w  constant 


9.  Tg 


(4-22) 
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Combining  equations  (4-21)  and  (4-22), 

Tad  To0’1  ,, 

— — I  «  constant  — (**-23) 

T9  P0* 

From  equation  (4-23)  it  is  seen  that  as  P  increases  the  gas  temperature 
Increases.  The  seemingly  contradictory  effects  of  increasing  heat  trans¬ 
fer  and  in  r  *asing  gas  temperature  result  from  the  coupling  of  the  heat 
transfer  rate  and  the  mass  flow  rate.  From  earlier  discussions  it  is  re¬ 
called  that  increasing  the  gas  temperature  produces  a  shift  in  the 
stability  limits  toward  lower  combustion  temperatures.  Hence  a  heat 
transfer  argument  leads  to  the  conclusion  that  increasing  P  produces  a 
shift  in  stability  limits  away  from  stoichiometric.  Since  this  conclu¬ 
sion  is  contrary  to  the  experimental  observations  it  is  felt  that  the 
pressure  dependence  of  heat  transfer  cannot  account  for  the  observed  shift 
in  stability  limits. 

Effect  of  P  on  Combustion  Zone  Analysis 

In  the  analysis  of  the  combustion  zone  boundary  condition 
pressure  and  temperature  were  related  through  the  isentropic  equations. 

As  noted  in  Chapter  3  the  isentropic  equations  are  not  valid  in  the 
combustion  zone  due  to  energy  release  and  diffusion.  In  employing  a 
pressure-dependent  energy  release  law  with  the  results  of  the  isentropic 
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coobustion  zone  analysis  Sirignano  and  Crocco  obtained  results  which  were 
not  in  agreement  with  experimental  observations .  To  circumvent  this 
problem  they  assumed  that  the  gas  rocket  combustion  process  could  be 
approximated  by  an  energy  release  law  which  was  independent  of  pressure 
(equation  3-1)*  It  is  not  surprising,  therefore,  that  the  theoretical 
stability  limits  are  independent  of  pressure. 

In  the  analysis  of  the  combustion  zone  boundary  condition 
Sirignano  and  Crocco  assumed  that  the  combustion  zone  thickness  (  &  ) 
was  constant.  As  noted  in  Section  C  J  is  a  function  of  time  and  de¬ 
pends  on  the  reaction  rate.  If  the  reaction  rate  were  pressure  dependent 
then  &  would  also  depend  on  pressure.  Thus  it  is  seen  that  to  properly 
account  for  pressure  effects  (for  an  Arrhenius -type  rate  function)  it  is 
necessary  to  include  non-isentropic  effects  and  a  variable  8  in  the 
combustion  zone  analysis. 

2.  Dependence  of  Instability  Strength  on  P 

To  determine  the  dependence  of  the  instability  strength  (A  P/P) 
on  the  mean  combustion  pressure  an  experiment  was  carried  out  in  which 
A  F/P  was  measured  as  a  function  of  $  (ie.  T)  for  several  different  P 
(for  given  L,  M  and  propellant  combination).  The  results  of  a  typical 
experiment  using  H  -air  are  presented  in  Figure  35.  In  this  figure  ^  P/P 
is  plotted  against  the  adiabatic  combustion  temperature  for  P  =  4.4  and 


11.2  atm.  The  data  were  obtained  for  L  ■  38  cm,  the  31-hole  shoverhead 
injector,  the  Cp  -  46  converging-diverging  nozzle  (M  =  0.015)  and  the 
oxidizer-rich  unstable  region.  The  heavy  vertical  bars  indicate  the 
variation  in  6hock  strength  during  a  given  run.  For  a  given  combustion 
temperature  the  instability  strength  for  P  «  4.4  atm  is  significantly 
less  than  the  instability  strength  for  P  ■  11.2  atm. 

In  the  Sirignano-Croceo  analysis  the  strength  of  the  instability 
is  independent  of  P  (due  to  the  fact  that  the  combustion  law  3-1  is  in¬ 
dependent  of  pressure).  In  Section  G-l  it  was  noted  that  a  pressure- 
independent  rate  law  was  necessary  due  to  certain  simplifications  in  the 
combustion  zone  analysis.  Ibe  data  of  Figure  35  are  presented  as  further 
evidence  for  the  importance  of  modifying  the  combustion  zone  analysis  so 
that  a  pressure-dependent  rate  law  can  be  employed. 

H.  Effect  of  Combustion  Chamber  Length  on  the  Instability 

I.  Dependence  of  Stability  Limits  on  L 

From  the  instability  maps  of  Figures  17  and  19  it  is  seen  that 

for  a  given  propellant  combination,  combustion  pressure,  injector  and 

exhaust  nozzle  the  value  of  can  depend  on  the  combustion  chamber 

limit 

length  (L).  Experimental  studies  of  longitudinal  combustion  instability 
in  gas  rockets  (27,  29,  31 )  have  shown  thet  for  a  given  propellant  com¬ 
bination  and  combustion  pressure  there  is  a  minimum  value  for  L  below 


which  the  combustion  is  stable.  Values  for  ^in  for  some  propellant  com¬ 
binations  of  interest  are  tabulated  below. 


;  Propellant 

Lmin  (cm) 

Source 

Hg-air 

12.7 

(27) 

Vair 

27.9 

(31) 

natural  gas -air 

18.0  i 

(34) 

In  the  Sirignano-Crocco  analysis  the  stability  limits  are  independent  of 
L  and  no  minimum  length  for  stable  combustion  is  predicted.  The  differences 
between  the  experimental  observations  and  the  theoretical  results  can  be 
explained  by  comparing  the  conditions  in  the  gas  rocket  with  the  assumptions 
of  the  theory.  From  the  discussion  of  Chapter  3  it  is  recalled  that  the 
Sirignano-Crocco  analysis  assumes  frictionless,  adiabatic  flow.  The  flow 
in  the  gas  rocket  is  highly  non-adiabatic  (See  Figure  1+3),  and  lrictional 
effects  are  present.  Earlier  it  was  shown  that  for  the  energy  release  law 
used  in  the  Sirignano-Crocco  analysis  losses  due  to  heat  transfer  and 
friction  tend  to  shift  the  stability  limits  toward  lower  combustion  temper¬ 
atures.  Heat  transfer  and  friction  losses  increase  with  increasing 
chamber  length.  Hence  as  L  increases  it  is  expected  that  the  stability 
limits  will  shift  toward  lower  combustion  temperatures- ie.  away  frot’. 
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stoichiooetric  .  This  shift  in  the  stability  limits  is  evident  for  the 

oxidirer-rich  unstable  regions  in  Figures  17  and  19.  The  variation  in 

the  adiabatic  corresponding  to  the  observed  change  in  $  limit 

for  H^-air  as  L  varies  from  19  cm  to  104  cm  is  tpproximatcly  130°K.  In 

the  vicinity  of  the  fuel- rich  unstable  region  for  H^-air  an  equivalent 

variation  in  T  corresponds  to  a  variation  in  <f>  of  0.04. 

limit  ^  limit 

This  variation  in  $  is  too  small  to  be  detected  experimentally, 

limit 

and  is  the  probable  explanation  for  the  apparent  lack  of  dependence  of 
$  nmit  on  ^  for  fuel-rich  unstable  regions. 

The  Sirignano-Crocco  analysis  assumes  that  the  combustion  zone 
thickness  is  significantly  smaller  than  the  chamber  length- ie.  that  the 
energy  is  added  to  the  oscillation  at  x  *  0,  In  the  gas  rocket  the 
combustion  zone  thickness  is  of  the  order  of  2-4  cm  (See  Section  K).  As 
the  chamber  length  decreases  below  40  cm  the  assumption  of  a  "zero¬ 
thickness"  combustion  zone  becomes  less  appropriate,  and  the  analysis 
becomes  subject  to  doubt.  In  addition  as  8  becomes  a  significant  frac¬ 
tion  of  L  energy  is  added  to  the  oscillation  in  a  non- optimum  manner.  Non¬ 
optimum  energy  addition  gives  a  steady-state  instability  strength  which  is 
less  than  the  steady-state  instability  strength  for  energy  addition  at  a 
pressure  antinode  (assuming  the  same  energy  release  rate  for  both  cases). 

As  the  chamber  length  decreases  the  combustion  zone  extends  over  a  larger 
fraction  of  the  chamber  and  hence  the  energy  added  to  the  oscillation 


decreases.  It  is  conceivable  that  for  sufficiently  short  chambers  (ie. 
for  pressure  oscillations  with  wavelengths  of  the  same  order  as  8  )  the 
energy  addition  would  not  be  sufficient  to  support  a  combustion  oscilla¬ 
tion.  This  conclusion  is  supported  by  the  minimum  length  data  of  Refer¬ 
ence  27  and  the  measurement  of  &  (Section  K)  for  H^-air: 


"min 


12.7  cm 


^  «  U  cm 


for  P  «  7-8  atm  and  a  10-hole  showerhead  'njector. 

2.  Dependence  of  Instability  Strength  on  L 

To  determine  the  dependence  of  the  instability  strength  (  A  P/P) 
on  the  length  of  the  combustion  chamber  (L)  an  experiment  was  carried  out 
in  which  A  P/P  waa  measured  as  a  function  of  $  (ie.  T)  for  several  values 
of  L  (for  a  given  P,  M  and  propellant  combination).  Ihe  results  of  a 
typical  experiment  using  H^-air  are  presented  in  Figure  36.  In  this 
figure^  P/P  is  plotted  against  the  adiabatic  combustion  temperature  for 
L  •  38  cm  and  76  cm.  The  data  were  obtained  for  P  »  U.U  atm,  the  31-hole 
6howerhead  injector,  the  rp  ■  U6  converging-diverging  nozzle  (M  «  0.015) 
and  the  oxidizer-rich  unstable  region.  The  heavy  vertical  bars  indicate 
a  non-periodic  variation  In  instability  strength  during  a  given  run.  For 


a  given  combustion  temperature  the  instability  strength  for  L  =  76  cm  is 
significantly  larger  than  for  L  ■  38  cm.  A  similar  length  effect  has  been 
reported  by  Zucrov  and  06born  ( 31 )  for  propane-air.  Some  results  from 
their  study  are  tabulated  below. 


propane-air 

<X>  =  1.2 

P  =  45  psia 


As  L  increases  from  12.6  inches  to  22.3  inches  A  P/P  increases  frotn  0.09 
to  O.U9. 


L  (in) 

A  p  (psi) 

12.6 

4 

16.3 

17 

22.3 

22 

28.0 

22 

In  Section  H-l  it  was  noted  that  due  to  non-optimum  i..ergy 
addition  the  energy  added  to  the  oscillation  decreases  as  the  chamber 
length  decreases.  It  follows  that  a6  the  chamber  length  decreases  the 
strength  of  the  oscillation  decreases.  Thus  it  is  seen  that  the  observed 
decrease  in  instability  strength  with  decreasing  L  can  be  explained  in 
terms  of  a  finite-size  combustion  zone,  and  is  dircetly  related  to  the 


existence  of  a  minimum  L  for  unstable  combustion. 
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An  alternative  explanation  for  the  length  dependence  of  the 
instability  strength  was  suggested  by  Zucrov  and  Osborn  (31).  Their 
explanation  is  based  on  an  assumed  increase  in  the  energy  available  to 
support  an  oscillation  with  increasing  L.  A  critical  analysis  of  the 
results  of  Zucrow  and  Osborn  by  Crocco  (32)  indicates  that  the  energy 
available  to  support  an  oscillation  is  relatively  insensitive  to  varia¬ 
tions  in  L.  Hence,  the  length  dependence  of  instability  strength  cannot 
be  explained  by  an  increase  in  available  energy  with  increasing  L. 

I.  Injector  Experiments 

1.  Efiect  of  Injector  on  Stability  Limits 

To  determine  the  effect  of  the  injection  pattern  on  the  stability 
limits  an  experiment  was  carried  out  in  which  the  regions  of  unstable 
combustion  for  I^-air  were  located  in  the  <£>  -L  plane  for  five  injectors 
with  different  injection  patterns.  The  injectors  used  in  the  experiment 
were:  porous  plug,  10-hole  showerhead,  31-hole  showerhead,  ^9-hole  shower- 
head  and  double-triplet  impinging.  A  detailed  description  of  these  in¬ 
jectors  is  presented  in  Chapter  2.  The  porous  plug  injector  provides  for 
nearly  one-dimensional  injection  of  propellants.  The  showerhead  injectors 
produce  injection  patterns  which  deviate  from  one-dimensional  injection. 

In  addition  the  combustion  zone  thicknesses  (  8  )  obtained  with  the 
showerhead  injectors  are  significantly  larger  than  the  8  obtained  with  the 


porous  plug  injector  (See  Section  K).  As  the  number  of  injection  orifices 
in  the  shoverhead  injectors  decreases  deviations  from  one-dimensional 
injection  and  the  combustion  zone  thickness  increase.  The  impinging 
injector  produces  a  strong  recirculation  pattern  and  hence  provides  an 
injection  pattern  which  deviates  significantly  from  one -dimensional  in¬ 
jection.  The  above  combination  of  injectors  permits  an  investigation  of 
the  effects  of  deviations  from  one-dimensional  injection  and  combustion 
zone  thickness  on  the  stability  limits. 

The  results  of  the  experiment  for  P  =  7-8  atm  and  the  p  =  46 
converging-diverging  nozzle  are  presented  in  Figures  30  and  37  through  39* 
From  Figures  30.  37  and  38  it  is  seen  that  as  the  number  of  injection 
orifices  in  the  shoverhead  injectors  increases  (ie.  as  the  injection 
pattern  approaches  one -dimensional  and  8  decreases)  the  stability  limits 
shift  toward  lower  combustion  temperatures.  For  the  porous  plug  injector 
stable  combustion  was  observed  over  the  entire  range  of  equivalence  ratios 
tested  (0.3  to  3-0).  This  result  is  consistent  with  the  data  of  Pelmas 
(27)  and  with  the  data  of  Figures  30»  37  and  38  (since  the  porous  plug 
injector  is  in  effect  a  "many-hole"  shoverhead  injector).  If  the  trend 
of  Figures  30,  37  and  38  is  extended  then  the  stability  limits  for  the 
"many-hole"  porous  plug  injector  should  fall  at  a  lower  combustion  temper¬ 
ature  than  the  stability  limits  for  the  49- hole  shoverhead  injector  (and 
thus  in  all  probability  outside  the  range  of  <X>  tested).  The  results  for 


the  impinging  injector  (Figure  39)  are  also  consistent  with  the  above 


observations.  The  impinging  injector  provides  a  non-uniform  injection 
pattern  and  a  diffuse  combustion  zone.  Hence  the  stability  limit  for 
the  impinging  injector  should  fall  at  higher  combustion  temperatures  than 
the  stability  limits  for  the  showerhead  injectors. 

The  important  observation  to  be  made  from  the  above  results  is 

that  the  dependence  of  the  stability  limits  on  the  injector  seems  contrary 

to  what  would  be  expected  from  the  Sirignano-Crocco  analysis.  The  porous 

plug  injection  closely  approximates  the  conditions  of  the  Sirignano-Crocco 

model  (ie.  one -dimensional  flow  and  zero- length  combustion  zone)  and  yet 

the  results  obtained  using  this  injector  are,  at  best,  in  poor  agreement 

with  the  theory.  If  the  stability  limits  fall  outside  the  range  of 

tested,  as  suggested,  then  the  value  for  E  would  differ  even  more 

limit 

from  the  accepted  values  for  the  activation  energy  than  is  the  case  for 
the  showerhead  injectors.  The  discrepancies  between  the  theoretical 
results  and  the  experimental  observations  probably  can  be  traced  to  the 
simplified  combustion  zone  analysis  employed  in  the  Sirignano-Crocco 
theory. 


2.  Effect  of  Injector  on  the  Location  of,  Regions  of  Harmonic 
Mode  Instabilities 


The  experiments  of  Section  1-1  also  provide  information  on  the 
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effect  of  the  injector  on  the  location  of  the  regions  of  harmonic  mode 
instabilities.  Of  particular  interest  are  the  experiments  using  the 
three  showerhead  injectors.  From  Figures  30,  37  and  38  it  is  seen  that 
for  given  <g>  the  minimum  critical  length  for  the  appearance  of  the  second 
harmonic  mode  decreases  as  the  number  of  injection  orifices  increases. 

This  result  can  be  explained  in  terms  of  the  harmonic  mode  criterion 
(4-19).  As  will  be  seen  in  Section  K  for  a  given  combustion  pressure  and 
equivalence  ratio  the  thickness  of  the  combustion  zone  decreases  as  the 
number  of  injection  orifices  increases-ie.  the  residence  time  of  a  gas 
particle  in  the  combustion  zone  (T)  decreases  as  the  number  of  injection 
orifices  increases.  From  (4-19)  the  lower  limit  for  a  given  harmonic 
mode  is  defined  by 


For  this  equation  to  be  satisfied,  as  T  decreases  0  must  decrease.  In 
other  words  as  the  number  of  injection  orifices  increases  the  minimum 
critical  length  for  a  given  harmonic  mode  decreases.  A  similar  argument 
can  be  advanced  to  explain  the  ...inimum  critical  length  for  the  second 
harmonic  mode  for  the  impinging  injector  (Figure  39)*  The  data  of  Section 
K  indicate  that  the  impinging  injector  gives  a  diffuse  combustion  zone. 
Hence  the  residence  time  (and  thus  0  )  increases.  The  argument  is  not 
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as  clear  cut  as  for  the  shoverhead  injectors  due  to  the  highly  non- 
uniform  injection  pattern  produced  by  the  impinging  injector. 

J.  Effect  of  Nozzle  Geometry  on  the  Instability 

1.  Dependence  of  Stability  Limits  on  Nozzle  Geometry 

To  determine  the  dependence  of  the  stability  limits  on  exhaust 
nozzle  geometry  the  regions  of  unstable  combustion  for  H^-air  were  located 
in  the  $  -L  plane  for  several  different  nozzles.  The  nozzleB  employed 
in  this  experiment  were:  two  converging-diverging  nozzles  with  convergence 
angles  of  60°  and  contraction  ratios  of  23  and  46  (Figure  12)  and  a  tung¬ 
sten  plug  nozzle  with  an  effective  contraction  ratio  of  46  (Figure  13). 

The  results  of  this  experiment  for  P  =  4.4  atm  and  the  31-hole  showerhead 
injector  are  shown  in  Figure  40.  Only  the  stability  limits  nearest  stoi¬ 
chiometric  for  the  oxidizer-rich  unstable  regions  are  shown. 

Hie  effect  of  the  length  of  the  convergent  section  of  the  nozzle 
on  the  stability  limits  is  obtained  by  comparing  the  limits  obtained  for 
the  smaller  converging-diverging  nozzle  (^p  ■  46)  and  the  plug  nozzle. 
These  two  nozzles  have  identical  contraction  ratios  (ie.  identical  flow 
areas)  and  convergent  sections  of  different  length  (The  plug  nozzle 
approximates  the  "zero-length"  nozzle  assumed  in  the  Sirignano-Crocco 
analysis).  From  Figure  40  it  is  seen  that  the  stability  limits  for  the 
two  nozzles  are  the  same  within  the  error  of  the  experiment  (Appendix  B). 
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In  view  of  the  results  of  Figure  40  it  is  concluded  that  the  assumption 
of  a  zero-length  nozzle  is  reasonable  for  nozzle  lengths  less  than  0.05  L. 

The  effect  of  nozzle  contraction  ratio  on  the  stability  limits 
is  obtained  by  comparing  the  limits  for  the  two  converging-diverging 
nozzles.  These  two  nozzles  have  identical  convergence  angles  and  differ¬ 
ent  contraction  ratios.  From  Figure  40  it  is  seen  that  the  stability 
limits  for  the  two  nozzles  are  significantly  different.  Ihe  stability 
limit  for  the  nozzle  with  the  larger  flow  area  (CR  =  23)  is  located  at 
lower  combustion  temperatures  than  is  the  limit  for  the  smaller  nozzle. 

I 

This  result  can  be  attributed  to  two  effects-  energy  losses  at  the  nozzle 
and  increased  heat  transfer.  The  energy  lost  frcm  a  wave  during  reflec¬ 
tion  increases  as  the  nozzle  flow  area  increases-ie.  the  stabilizing 
effect  of  the  nozzle  increases  as  the  flow  area  increases.  Hence  as  the 
contraction  ratio  decreases  the  stability  limits  should  shift  toward 
lower  combustion  temperatures.  To  maintain  a  given  P  the  mean  mass  flow 
rate  increases  as  the  nozzle  flow  area  increases.  As  noted  earlier  in¬ 
creasing  the  mass  flow  rate  increases  the  heat  transfer  to  the  combustion 
chamber  walls.  Hence,  as  before,  the  stability  limits  should  shift 
toward  lower  combustion  temperatures. 

2.  Effect  of  Nozzle  Geometry  on  Instability  Waveform 


To  determine  the  effect  of  the  exhaust  nozzle  on  the  instability 
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waveform  pressure-time  histories  at  a  station  Just  upstream  from  the 
nozzle  entrance  plane  were  obtained  for  three  different  nozzles-  the 
tungsten  plug  nozzle,  the  Cp  =  46  converging-diverging  nozzle  and  the 

n 

p  ■  23  converging-diverging  nozzle.  Four  typical  pressure-time  histor¬ 
ies  are  presented  in  Figures  4l  and  42.  These  histories  were  obtained 
using  H2-air  for  L  =  104  cm,  P  *  6.1  atm,  critical  nozzle  flow,  the  31- 
hole  shoverhead  injector  and  the  oxidizer-rich  unstable  region.  The 
pressure  transducer  was  flush-mounted  in  the  wall  of  the  combustion  cham¬ 
ber  approximately  5  cm  upstream  from  the  nozzle  entrance  plane.  The  mode 
of  the  instability  is  the  second  harmonic.  Figure  4la  is  a  typical 
pressure-time  history  for  the  tungsten  plug  nozzle.  The  first  pressure 
rise  corresponds  to  a  wave  incident  on  the  nozzle;  the  second  pressure 
rise  corresponds  to  the  wave  after  reflection  from  the  nozzle.  It  is  seen 
that  both  the  incident  and  reflected  waves  are  well-defined  shock  waves 
and  that  no  significant  loss  in  shock  strength  occurs  upon  reflection. 
Figure  4lb  is  a  typical  pressure-time  history  for  the  CR  =-  46  converging- 
diverging  nozzle.  This  nozzle  has  the  same  flow  area  as  the  plug  nozzle. 
It  is  seen  that  some  distortion  of  the  wave  occurs  on  reflection  from  the 
nozzle  but  that  the  reflected  wave  remains  shock-like.  No  significant 
loss  in  shock  strength  occurs  upon  reflection.  At  a  given  combustion 
temperature  there  is  no  significant  difference  between  the  shock  strengths 
for  the  plug  nozzle  and  the  Cp  =  46  converging-diverging  nozzle.  Figure 
42a  is  a  typical  pressure-time  history  for  the  p  =  23  converging- 


diverging  nozzle.  This  nozzle  has  the  seune  convergence  angle  as  the  ^  = 
46  converging-diverging  nozzle  and  twice  the  flow  area.  It  Is  seen  that 
the  reflected  wave  Is  distorted,  but  still  shock- like,  and  that  there  is 
a  measurable  loss  in  shock  strength  on  reflection.  From  the  results  of 
this  experiment  it  is  concluded  that  for  a  given  oscillation  frequency 
(ie.  L)  the  loss  in  shock  strength  on  reflection  depends  primarily  on 
nozzle  flow  area  and  is  relatively  independent  of  the  length  of  the  con¬ 
vergent  section  of  the  nozzle.  A  similar  result  was  obtained  by  Feller 
(47)  in  a  shock  tube  study  of  the  reflection  of  shock  waves  from  critical 
flow  nozzles.  Although  some  distortion  of  the  reflected  shock  front  1b 
observed  with  the  converging-diverging  nozzles  the  wave  remains  shock-llke 
and,  in  fact,  soon  steepens  into  a  sharp  shock  front.  Figure  42b  is  a 
typical  pressure-time  history  at  a  station  15  cm  upstream  from  the  en¬ 
trance  plane  of  the  cp  =  46  nozzle.  Note  that  the  distorted  wave  front 
of  Figure  4lb  has  steepend  into  a  well-defined  shock  front. 

K.  Measurement  of  Axial  Temperature  Profiles 

In  several  of  the  experiments  described  in  this  chapter  axial 
distributions  of  mean  gas  temperature  were  obtained  using  a  series  of 
thermocouples  located  on  the  centerline  of  the  combustion  chamber.  A 
detailed  description  of  the  thermocouples  is  presented  in  Chapter  2. 

The  temperature  distributions  are  used  to  determine  the  mean  combustion 


temperatures  at  the  stability  limits  (T  )  and  the  thickness  of  the 

limit 

combustion  zone  (  g  ).  In  addition  the  temperature  distributions  give  a 
quantitative  measure  cf  the  total  heat  transfer  from  the  combustion  gases 
(See  Assumption  5>  Chapter  3)» 

In  Figure  43  typical  temperature  profiles  are  shown  for  the 
various  injectors  used  in  the  experiments.  The  profiles  shown  were 
obtained  using  H^-air  with  7  -  7*8  atm,  L  *  104  cm  and  the  plug  nozzle. 

The  equivalence  ratio  and  corresponding  adiabatic  combustion  temperature 
are  given  for  each  profile. 

Before  considering  the  results  of  Figure  43  several  details  of 
the  measurement  procedure  will  be  discussed.  The  thermocouples  used  in 
the  experiments  were  not  equipped  with  radiation  shields.  To  obtain  an 
estimate  of  radiation  losses  an  experiment  was  carried  out  in  which  the 
gas  temperature  at  a  given  location  in  the  combustion  chamber  was  measured 
simultaneously  by  a  radiation  shielded  thermocouple  and  an  unshielded 
thermocouple.  The  results  of  this  experiment  indicated  that  within  the 
scatter  of  the  temperature  data  (il5°K)  there  was  no  difference  between 
the  two  thermocouple  readings.  This  result  is  in  agreement  with  the  re¬ 
sults  of  a  heat  transfer  analysis  of  the  thermocouple  system.  The  fact 
that  the  readings  of  the  shielded  and  unshielded  thermocouples  are  approx¬ 
imately  the  same  is  due  to  the  very  large  convective  heat  transfer 


coefficient  of  the  combustion  gases.  Since  it  appears  that  no  correction 
for  radiation  losses  is  necessary  the  temperature  data  of  Figure  43  are 
direct  thermocouple  readings.  Temperature  measurements  were  made  using 
uncoated  platinum- platinum/rhodium  thermocouples.  It  is  well  known  that 
platinum  lias  a  catalytic  effect  on  the  hydrogen- oxygen  reaction.  Hence 
the  temperature  measurements  made  in  the  reaction  zone  itself  cure  probably 
not  reliable.  Downstream  from  the  reaction  zone  (ie.  downstream  from  the 
peak  temperature)  the  temperature  data  should  be  reliable  to  within  ii£ 
(the  accuracy  of  the  thermocouple).  The  thermocouple  configuration  is 
such  that  the  temperature  measured  is  a  stagnation  temperature.  Since 
the  mean  Mach  number  of  the  combustion  gases  is  very  small  (  <  0.03)  the 
stagnation  and  static  temperatures  are  for  all  intents  and  purposes  the 
same. 


From  Figure  43  it  is  seen  that  for  a  given  equivalence  ratio 
the  peak  temperature  is  essentially  independent  of  the  injector  and 
approximately  40°K  less  than  the  adiabatic  combustion  temperature.  This 
result  is  not  surprising  since  the  heat  losses  from  the  relatively  thin 
combustion  zone  are  small.  It  is  seen  that  the  combust:  n  zone  thickness 
for  the  porous  plug  injector  is  significantly  less  than  the  8’s  for 
the  showerhead  and  impinging  injectors.  Furthermore  8  decreases  as  the 
number  of  injection  orifices  in  the  showerhead  injectors  increases. 
Finally  it  is  noted  that  the  decrease  in  gas  temperature  over  the  length 
of  the  combustion  chamber  (*wlOO  cm)  is  considerable  ( ~  700°K)  so  that 
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the  aa sumption  of  adiabatic  flow  (for  lorg  motors  at  least)  is  not 
particularly  good. 
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CHAPTER  5.  CONCLUSIONS 


A.  Mechanism  for  High-Frequency  Longitudinal  Gas -Phase 
Combust ion  Instability 

The  Sirlgnano-Croccc  combustion  instability  analysis  is  based 
on  a  gas-phase  chemical  kinetic  driving  mechanism,  and  employs  a  simpli¬ 
fied  model  for  the  gas  rocket  combustor  (Chapter  3)  •  In  comparing  the 
experimental  observations  vlth  the  theoretical  results  (Chapter  4)  it  was 
found  that  within  the  limitations  Imposed  by  the  analytical  assumptions 
there  is  relatively  good  agreement  between  theory  and  experiment.  The 
theory  provides  a  criterion  for  the  boundaries  separating  regions  of 
stable  and  unstable  combustion  (equations  3-4  and  4-6)  vhich  was  verified 
qualitatively  by  experiment.  In  addition  the  theory  gives  the  correct 
dependence  of  instability  strength  on  the  mean  Mach  number  of  the  combus¬ 
tion  gases  and  the  mean  combustion  temperature.  The  several  areas  of 
disagreement  between  theory  and  experiment  generally  can  be  explained  in 
terms  of  differences  between  the  conditions  in  the  gas  rocket  and  the 
theoretical  model.  The  results  of  Chapter  4  indicate  that  the  most  im¬ 
portant  differences  between  theory  and  experiment  can  be  traced  to  over¬ 
simplifications  in  the  combustion  zone  boundary  contrition  and  the  energy 
release  lav.  The  experimental  results  point  up  the  fact  that  the  effects 
of  a  finite  combustion  zone  on  the  instability  phenomenon  are  not  fully 
understood.  Crocco  has  recently  advanced  an  instability  theory  (33) 
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vhlch  Incorporates  a  more  realistic  energy  release  lav  and  a  combustion 
zone  of  finite  thickness.  Hardesty  (49)  is  presently  carrying  out 
experiments  to  better  understand  the  effect  of  a  finite  combustion  zone 
cn  the  instability  phenomenon. 

Various  experimental  studies  of  longitudinal  combustion  in¬ 
stability  in  gas  rockets  have  been  carried  out  (27-29,  31,  34).  In 
Chapter  4  and  Reference  50  the  results  of  these  studies  have  been  shown 
to  be  consistent  with  the  results  of  the  present  investigation.  Avail¬ 
able  gas  rocket  data,  together  with  the  results  of  the  Sirignano-Crocco 
analysis,  provide  strong  evidence  for  a  mechanism  for  high-frequency 
longitudinal  gas-phase  combustion  instability  based  on  the  interaction 
between  the  combustion  oscillation  and  the  energy  release  rate-ie. 
chemical  kinetics. 

B.  Possible  Applications  of  Gas  Rocket  Instability  Data 

The  most  obvious  application  of  gas  rocket  instability  data  is 
to  homogenous  gas-phase  chemical  kinetics.  If  a  theory  could  be  developed 
which  provides  a  quantitatively  accurate  criterion  for  the  stability 
limits  in  terms  of  the  chemical  kinetic  parameters  then  it  would  be 
possible  to  use  stability  limit  data  obtained  from  gas  rocket  experiments 
to  determine  the  overall  kinetics  for  the  gas -phase  reaction  of  a  given 


propellant  combination 
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If  a  comparison  is  made  between  the  experimentally  observed 
properties  of  the  longitudinal  mode  of  high-frequency  combustion  insta¬ 
bility  in  gaseous  propellant  and  liquid  propellant  motors  it  is  found 
that  in  a  number  of  ways  they  are  similar.  Under  certain  conditions 
the  instability  waveform  are  identical-ie.  shock  discontinuities  followed 
by  exponential  decay  in  pressure  and  velocity,  with  oscillation  frequencies 
corresponding  to  one  of  the  resonant  frequencies  of  the  combustion  chamber. 
The  nature  and  ordering  of  the  regions  of  harmonic  mode  instabilities  are 
the  same.  In  spite  of  these  similarities  it  is  felt  that  due  to  several 
uncertainties  in  the  gas-phase  instability  phenomenon  an  attempt  to  re¬ 
late  gas-phase  and  two-phase  combust i»n  instaoility  would  be  premature. 
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APPENDIX  A:  CALIBRATION  PROCEDURES 


Critical  Flow  Orifices 


The  mass  flow  rate  of  a  perfect  gas  through  a  critical  flow 
orifice  is  given  by 


m  =  g(r,*) 


(a-i) 


where  m=  mass  flow  rate,  g(  y  function  of  the  specific  heat  ratio 

and  molecular  weight,  Aorifice=  exit  area  of  the  orifice,  P^=  stagnation 

pressure  and  T  =  stagnation  temperature.  For  the  flow  of  a  real  gas 
o 

through  a  critical  flow  orifice  equation  (A-l)  must  be  modified  to  in¬ 
clude  the  effects  of  friction  and  compressibility-  ie. 

rns^l-flOr.’,,)A0),|F|CE-^-  (A-2 

where  C^=  discnarge  coefficient  <  1  (a  function  of  the  Reynolds  number) 

and  y,  ~  compressibility  factor.  It  is  noted  that  in  equation  (A-l) 

Y  is  a  constant  (for  a  given  gas)  whereas  in  equation  (A-2)  ^  is  a 

function  of  f  .  The  compressibility  factor  and  Y  are  tabulated  as 

functions  ol  P  for  several  gases  in  Reference  48.  The  discharge 
o 

coefficient  is  not  known  a  priori  so  that  it  is  not  possible  to  obtain 
solutions  to  equation  (A-2).  To  determine  m  as  a  function  of  P  and  To 
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the  critical  flow  orifices  used  in  the  gas  rocket  propellant  feed  system 
were  calibrated  using  a  water  displacement  technique.  In  this  technique 
the  mass  flow  rate  (for  known  and  Tq)  is  determined  by  allowing  the 
gas  to  displace  water  from  a  container  (see  sketch).  By  measuring  the 
volume  of  the  water  displaced  during  a  known  time  and  the  pressure  and 
temperature  of  the  gas  in  the  container  it  is  possible  to  calculate  the 
mass  of  gas  which  flowed  into  the  container  during  the  known  time  (ie. 
the  mass  flow  rate) . 


PRESSURE  GAUGES 


SKETCH  OF  ORIFICE  CALIBRATION  RIG 
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By  performing  this  measurement  for  several  Pq  and  referencing  the  result 
to  a  fixed  stagnation  temperature  (T  raf)  a  calibration  curve  for  the 
orifice  for  a  given  gas  and  T0  ref  is  obtained.  This  calibration  curve 
may  be  extended  to  other  Tq  by  using 


rhCl^)*m(T0R£p)  (a-3) 

T0  REF 

Nine  critical  flow  orifices  were  calibrated  using  Ng  and/or  He  in  the 
pressure  range  of  0-1800  psig  (1-123  atm-abs).  The  resulting  calibration 
curves  were  extended  to  the  propellant  gases  using 


m 


PROP 


m 

CAL 


(A-4) 


with  tabulated  values  for  Y  and  5J,  (Reference  48)  and  experimental 
values  for  Cp.  A  discussion  of  experimental  uncertainties  in  the 
calibration  data  is  presented  in  Appendix  B. 


Pressure  Transducers 


At  regular  intervals  the  steady-state  and  transient  pressure 
transducers  were  calibrated  using  a  dead  weight  testor.  The  output  from 
a  given  stead} -stute  transducer  was  recorded  on  the  chart  recorder  used 
to  record  the  transducer  output  during  the  experiments  so  that  the 


complete  pressure  recording  system  was  calibrated.  The  transient  pressure 
transducers  were  calibrated  statically,  recording  the  transducer  output 


on  the  visicorder. 


APPENDIX  B:  EXPERIMENTAL  UNCERTAINTIES 


edifice  Calibration  Data 

The  major  uncertainty  in  the  orifice  calibration  data  is 
introduced  in  extending  the  data  from  the  calibration  gases  (N,,,  He)  to 
the  propellant  gases  (H2>  CH^,  CO,  0o,  air).  This  extension  requires 
a  knowledge  of  the  ratio  of  the  discharge  coefficients  for  the  propellant 
gas  and  calibration  gas  (see  equation  A-4) .  In  using  equation  (A-U)  it 
was  assumed  that  the  Reynolds  number  dependence  of  C^  was  the  same  for 
both  the  propellant  gas  and  the  calibration  gas  (The  Reynolds  number 
dependence  of  C^  for  the  calibration  gas  was  obtained  from  the  calibra¬ 
tion  data).  The  uncertainty  in  the  calibration  data  is  obtained  by 
combining  estimates  of  the  uncertainty  introduced  by  the  C^  assumption 
with  the  uncertainty  in  the  measured  parameters  of  the  calibration 
experiment.  The  estimated  uncertainties  for  several  sets  of  calibration 
data  are  tabulated  below. 


orifice 

propellant 

uncertainty  ($) 

7 

H2 

±1.5  ! 

6 

CO 

+0.5 

3 

air 

+0.5 
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Equivalence  and  Dilution  Ratios 

Two  factors  contribute  to  the  uncertainty  in  the  equivalence 
and  dilution  ratios.  The  first  factor  is  the  uncertainty  in  the  orifice 
calibration  data.  The  second  factor  is  the  uncertainty  in  the  measure¬ 
ment  of  the  experimental  variables  (stagnation  pressure  and  temperature). 
Combining  these  two  factors  the  uncertainty  in  the  equivalence  and 
dilution  ratios  for  the  various  propellant  combinations  were  estimated. 
The  results  are  tabulated  below. 


propellant 

combination 

uncertainty  in  $  ($) 

uncertainty  in 
dilution  ratio  ($>) 

4>  >  i 

<J>  <  1 

I^-air 

*2.5 

*1.5 

I^-dilute  air 

*2.5 

-1-5 

*1.0 

CHj^-air 

-2.0 

-1*5 

*1.0 

CHu-(.U02-.6N2) 

*2.0 

-1.5 

*1.0 

CO-air 

-1*5 

-1.5 

• 
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It  is  noted  that  uncertainties  in  the  equivalence  and  dilution 
ratios  result  in  corresponding  uncertainties  in  the  adiabatic  combustion 
temperature  (eg.  for  H^-air,  $  =  2.50±0.06  gives  Tad=  lU25±25  °K). 

Transient  Pressure  Data 

The  frequency  response  of  the  transient  pi-essure  recording 
system  was  flat  to  within  10$  from  20  cps  to  6  kc  and  to  within  20$  from 
15  cps  to  12  kc.  To  determine  the  effect  of  the  inherent  filtering  of 
the  recording  system  on  the  experimental  data  a  1  kc  square  wave  was 
fed  into  the  recording  system.  The  output  and  input  are  shown  below. 


It  is  seen  that  there  is  only  a  very  slight  difference  between 


the  input  waveform  and  the  output  waveform.  It  is  concluded,  therefore, 


that  for  oscillation  frequencies  less  than  1  kc  the  effect  of  inherent 


filtering  on  the  transient  pressure  data  is  negligible. 


Temperature  Measurement 


See  Chapter  U. 


C-l 


APPENDIX  C:  SOME  RESULTS  OF  THE  SIRIQWABO-CROCCO  THEORY 

The  strength  of  the  shock  wave  instability,  in  terms  of  the 
changes  in  gas  velocity  and  pressure  across  the  wave,  is  given  by  (23) 


tests-™'’ 


* 

f 


AuBC_€Z  I- exp  (-X) 

~ x  TZtvTlJ  +  ff{€ z) 


(C-2) 


l'«KP 

l+exp 


(c-3) 


AP 

x  l+exp  f  X)  *  ° ' 


(C-4) 


It  is  noted  that  the  strength  of  the  instability  is  independent  of 
time  sind  is  altered  only  by  reflection  from  the  nozzle  or  injector. 
As  H  -•»  0,  */  0  and  thus  the  changes  in  u  and  P  across  the  wave 

are  the  same  for  both  the  downstream  and  upstream  moving  waves. 


The  velocity  of  the  shock  wave  is  given  by  (23) 


L 


C-2 


^.o+R)+fi 


(c-5) 


(c-6) 


It  is  noted  that  the  shock  velocity  increases  in  time  and  is  altered 
by  reflection. 


The  period  of  the  oscillation  is  given  by  (23) 


te)Ua£yii)!}  • 


{f-s$  -  ■  }*«* 


(C-7) 


It  is  noted  that  the  period  of  the  oscillation  corresponds  to  the 
fundamental  longitudinal  acoustic  mode  of  the  chamber  with  a  correction 
for  a  finite  strength  wave  and  mean  flow. 


The  equations  for  the  gas  velocity  and  pressure  behind  the 
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downstream  moving  shock  (as  functions  of  X  and  t  )  are  (23) 

?•■•?{  -i'/Jff?'1  •  [»"[»■•<».]  -  <«. 

«P  -  Tv  }  +  *&) 

exp  [--£(UM)x]j-  j-  ♦  &  (e2) 


It  is  noted  that  both  prescure  and  velocity  deca^  exponentially  in  time 
behind  the  shock  vave. 


The  equations  for  the  gas  velocity  and  pressure  behind  the 
upstream  moving  shock  (as  functions  of  X  and  t  )  are  (23) 


' “+ -X"  {tXP|^)>3  ~[(pf)  •*[*<■-«>»]- 

expXexp[-|(|+ra))(]j-jlL|+  O'Ur2) 


(c-10) 
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=  1  + 


^)«xp[-2  (•-«)*]  + 
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APPENDIX  D:  GAS  DATA 


Gas  I  Source  I  Grade 


General 

Hg  Dynamics  Electrolytic 


CF^  I  Matheson  I  Conmercial 


CO  I  Matheson  I  C.  P 


0  Air  Products  Dry 

2  ; 


Air  Reduction 


Liquid 


Typical 

Analysis 

Gas 

Concentration 

H 

2 

99-3*  (vol> ) 

°2 

0.2* 

Dev  point 

not  available 

93.6*  (mole) 

C2% 

3-6* 

% 

1.0* 

co2 

0.7* 

N2 

0.5* 

CO 

99-5*  (vol.) 

C°2 

200  ppm 

°2 

20  ppm 

N2 

75  PP«n 

HgO 

80-130  ppm 

°2 

99.5*  (mole) 

*a 

0.05* 

Ar 

0.45* 

Dew  point 

-76°F 

n2 

99.95*  (vol.) 

^  0 

Dev  point 

-80  F 
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